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Abstract 

These notes are a fnller version of fonr lectnres given at the 2015 Inter¬ 
national Snmnier Workshop in Reaction Theory held at Indiana University, 
Bloomington. The aim is to provide a simple introdnction to how the tools 
of “the A-matrix era” - i.e. the constraints of nnitarity, analyticity and cross¬ 
ing symmetry - can be incorporated into analyses of hnal state interactions 
in two- and three-hadron systems. The main focns is on corrections to the 
isobar model in three-hadron hnal states, which may be relevant once more 
as much larger data sets become available. 
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1 Introduction 


These lectures aim to give a simple introduction to the application of uni- 
tarity, analyticity and crossing symmetry - the main principles of S' - matrix 
theory (Eden et al. [1]) - to the analysis of hnal state interactions in two- 
and three-hadron systems. 

S - matrix theory flourished in the late 1950s and on through the 1960s. 
It was developed as a theory of the strong interactions between hadrons, 
to which the perturbative procedures of quantum held theory seemed inap¬ 
plicable. It is fair to say that S - matrix theory had only limited success 
as a hrst principles technique for calculating strong interaction amplitudes, 
though some important products survive, such as Regge theory (and re¬ 
markably enough S - matrix theory gave birth to string theory). Of course, 
strong interactions came to mean QCD, where both perturbative and non- 
perturbative (lattice) techniques have been very successful. Nevertheless, 
ab initio calculations of few hadron dynamics present a challenge, and S- 
matrix principles remain as valid constraints which should be incorporated 
into phenomenological analyses. 

Although some features, such as isolated resonances, show up clearly 
on simple intensity plots, in many cases we are interested in more subtle 
questions related to phases of amplitudes. Such information will have to 
come, as usual in quantum mechanics, from interferences. I briefly outline 
two (oversimplified) examples. 

Suppose we want to study two excited nucleon states and N 2 which 
are produced from an initial ttN state, and which decay sequentially to ttttN 
via the two decay chains —)■ ttA — )■ ttttN and N 2 —)■ pN —)■ ttttN. Then a 
simple model (essentially the isobar model) for the amplitude leading to the 
hnal ttttN state takes the form 

F = Aiti -\- A2t2 (1) 

where the Ai are the strong production amplitudes for N)| and N 2 , and the ti 
are the two-body hnal state interaction amplitudes in the A and p channels. 
Then |Fp will contain an interference term proportional to cos(0-|-0), where 
9 is the relative phase of Ai and A 2 , and (j) is the relative phase of ti and 
t 2 - So from the intensity |Fp we can learn about the relative phase of the 
production amplitudes, provided that we know the relative phase of the hnal 
state two-body amplitudes ti and ^ 2 - In the isobar model, these are assumed 
to be determined from the known two-body scattering data. But we will see 
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that unitarity (or equivalently rescattering amongst the hnal state particles) 
forces corrections to the isobar model, which affect these relative phases. 
At some point, therefore, such corrections should be incorporated into the 
analysis. 

A second example concerns the extraction of CP-violating phases in states 
decaying weakly to hadronic hnal states. Again we can use (1) to make the 
point, where now Aj are the weak production amplitudes and as before the 
U are strong two-body hnal state amplitudes. The CP-conjugate amplitude 
will be 

F = AlF + Alt2, (2) 

and the CP-violation will be observable from the diherence 

|F|2 - |F|2 = 4 Im(AiA;) lm{t2tl). (3) 

To get an ehect, there needs to be a phase diherence between both the two 
weak amplitudes and the two strong amplitudes. And to extract the value 
of the CP-violating weak phase diherence we need to be sure of the strong 
phase diherence. In two-body hnal states the latter is known from two-body 
data, but in three-body states rescattering ehects will again modify the U 
phases. 

It would be nice if we could have a phenomenology that was independent 
of approximations necessarily made in describing the hadronic hnal state 
interactions. Such a model-independent analysis generally requires very large 
data samples. Although these may now be beginning to be available, it seems 
likely that amplitudes with some theory ingredients will still be needed for 
some time. And with vastly more data, the dehciencies in models like the 
isobar model may need to be remedied. A reasonable way to tackle this is to 
require as a “minimum theory” that our amplitudes satisfy the old S - matrix 
principles mentioned previously - that is, we aim to provide amplitudes which 
at least obey the constraints of unitarity, analyticity, and crossing symmetry, 
as far as possible. 

These lectures will describe what these constraints are and how they 
are implemented in some simple examples. We begin with two-hadron hnal 
states, introducing unitarity and the iP-matrix. Then we add analyticity, and 
dispersion relations. Our main focus, though, will be on three-hadron hnal 
states. We show how unitarity in the two-body sub-energy channels places 
a constraint on the isobar decay amplitudes, and how analyticity enables 
us to convert this into integral equations for modihed isobar amplitudes. 
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which satisfy two-body unitarity. We shall see that, somewhat surprisingly, 
these amplitudes can actually satisfy three-body unitarity as well. This “two- 
body” approach to what is after all a three-body problem is conceptually very 
simple, and produces amplitudes which can directly replace the conventional 
isobar amplitudes. The price to be paid en route is a certain amount of 
gymnastics in the complex plane. 

Throughout we shall restrict ourselves to the simplest possible spin and 
angular momentum conhgurations, so that the logic of “unitarity -|- analyt- 
icity + crossing symmetry” can be clearly exhibited, unencumbered by other 
complications. However, I shall briefly report on the results of calculations 
from the 1970s and 1980s made for various physically realistic three-hadron 
systems. But this is not a review: rather, it mostly describes work that I 
was myself involved with, and no attempt is made to be comprehensive. 

2 Elastic 2—> 2 Unitarity 

2.1 One channel, one resonance 

2.1.1 Unitarity 

The unitarity relation for the T-matrix is 

T-T^ = 2iTpT^ = 2iT^pT (4) 

where p is the appropriate intermediate state phase space. For simplicity we 
consider the elastic scattering of two identical spinless bosons of unit mass, 
interacting in the I = 0 partial wave only. Then (4) becomes 

T(i,)-r(s) = 2v(i,)|r(s)p (5) 

or equivalently 

Imr(s) = p(s)|r(s)p (6) 

where s = 4 -|- 4g^ is the square of the total c.m. energy, q is the c.m. 
momentum, and (in a convenient normalization) 

pW=(^)“^ (7) 
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Figure 1: Resonance amplitude 


More generally, for a two-body threshold with unequal masses mi and m 2 
the phase space would be 


1 

s 


{mi + 7712 )^][s - 


(mi - m2) 


211/2 


-k{s, ml, 
s 



( 8 ) 


where 

k{a, b, c) = — 2ab — 2bc — 2ca. (9) 

A parametrisation satisfying (6) is 


r(s) = e“^sinh/p = (pcot(5 — ip) ^ 


( 10 ) 


where S is the phase shift. In particular, if we choose 6 = tan“^[pp^/(s^ — s)] 
the phase shift will rise from zero at threshold to vr as s ^ cxo, passing through 
7r/2 at s = Sr- This is a standard Breit-Wigner type resonance formula, with 
amplitude 


/(s) 



( 11 ) 


Near the peak of a narrow resonance we may set p(s) ~ p(sr); then the 
resonance maximum is reached at s = Sr, and the full width at half height 
in a plot of |T(s)p versus s is 2p(sr)p^. We represent f{s) by hgure 1. 


2.1.2 The complex plane, Reimann sheets 

We shall frequently be considering variables such as s to be complex, and 
our amplitudes will be assumed to be analytic functions of their arguments^. 
An example immediately arises in the case of the function p(s) in (4). In 

^Burkhardt’s book [2] contains a useful long first section on complex variable analy¬ 
sis, and then continues with equally relevant sections on collision theory and S'-matrix 
dynamics. 
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arg(s — 4) 


s = 0 



Figure 2: The branch points and cuts of p{s) 


the unitarity relation as written in (4) it is implicit that s > 4, the elastic 
scattering threshold. In that case, p{s) should be multiplied by 6*(s —4). But 
this is not an analytic function of s. Rather, we shall understand (4) to be 
true as it stands, and allow p{s) to be dehned for all values in the complex 
s plane, by analytic continuation from the physical region. That region is 
the real axis s > 4, approached from above: lim^^o /(s + ie). We need to be 
careful how we approach the real axis because, as we now discuss, it makes 
a difference due to the singularity structure of p{s). 

Viewed as an analytic function of the complex variable s, p{s) has branch 
points at s = 0 and s = 4, with associated cuts as shown in figure 2. The 
physical region for our 2-particle scattering is the real axis s > 4. Suppose 
we start at a point just above the real axis, with s > 4, with the square 
root function defined to be positive. We can continue the function \/s — 4 
on a circular path encircling the point s = 4, starting at a point just above 
the real axis s > 4, passing between s = 4 and s = 0, and returning to the 
real axis to the right of s = 4 but just below the real axis. On the real axis 
in the region 0 < s < 4 the square root becomes i\/4 — s, and at the end 
of the trip it has become —\/s — 4. Notice that the value of the function 
\/s — 4 for s > 4 and just above the real axis is not the same as the value 
of y/s — 4 for s > 4 and just below the real axis. That is why we draw a 
“cut” along the real axis s > 4, to remind ourselves of this discontinuity in 
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the function a/(s — 4). It is the value reached from above the real axis that 
is the “physical limit”. 

In addition to the branch point at s = 4, p{s) also has a branch point at 
at s = 0. Whereas the branch point at s = 4 has a clear physical origin - 
namely the two-particle threshold - that at s = 0 does not: it may be called 
a “kinematic” singularity. We shall see in section 4.1 how to get rid of it. 
We therefore continue to focus on the square root branch point at s = 4. 

We have discussed making one complete circuit of the point s = 4, ending 
up just below the real axis, to the right of s = 4. Let’s continue on from 
this point, and make another complete circuit. On this second circuit, the 
argument of s — 4 starts at the value 27r, and ends at the value 47r. Half-way 
round this second circuit, the argument has the value Stt. For the square 
root, we have to halve the argument, so the square root function starts at 
the value —\/s — 4, becomes —i\/4 — s half way round, and ends at the value 
y/s — 4 after the complete (second) circuit. Thus after two complete circuits 
around s = 4 the square root function returns to its original value. 

This description has been in terms of a double-valued function ii/s — 4 
dehned over a single complex plane. The standard alternative description 
shifts the multi-valuedness from the function to the space over which it is 
dehned. In the present case, we will have two complex planes, called “sheets”. 
On the hrst sheet, we use the positive square root -Ti/s — 4, and on the second 
sheet we use the negative square root —y/s — 4. On each sheet, we are dealing 
with a single-valued function. The interesting thing is that the sheets are 
connected in the region of the cut. Going once around s = 4 on sheet I, say, 
the function y/s — 4 ends up at the value —y/s — 4, which is just the same as 
the value of the function on the second sheet (namely —y/s — 4) evaluated 
just above the cut. So after one revolution on sheet I we pass smoothly onto 
sheet II as we cross the real axis. Continuing round on sheet II, we arrive 
after one circuit at a point just below the real axis s > 4, where the second 
sheet function takes the value -|-\/s — 4. This is the same as the value we 
started with in sheet I, before the two circuits. So the second time we cross 
the axis to the right of s = 4, we are back on sheet I. The way the sheets are 
connected along the cut is indicated in hgure 3. 

Another description of the square root function is also possible, and per¬ 
haps easier to visualize. Since the “sheet” business is all to do with the 
square root, maybe things would be simpler if we introduced a new variable 
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Figure 3: The two-sheeted function i/s — 4: a circuit on sheet I followed by 
a circuit on sheet II gets back to sheet 1. 

which is the square root itself, rather than s: namely, dehne a new variable 



( 12 ) 


which is of course just the magnitude of the momentum. If we set s = 
4 -|- re‘®, which parametrises a circle with centre at s = 4 and radius r, then 
q = It follows that the whole of the first s-sheet with 0 < 6* < 27r 

corresponds to the upper half g-plane, while the whole of the second s-sheet 
(27r < 6 < dvr) corresponds to the lower half g-plane. Two trips around the 
s = 4 threshold in s correspond to just one trip around g = 0, as shown in 
hgure 4. 

2.1.3 Resonance poles 

There is nothing more to be said about the square root function. What 
about the singularities of /(s) as given by (11)? It will simplify matters if 
the branch point at s = 0 is not present. We will see how to get rid of it 
in section 4.1, but for the moment it will be sufficient to replace the phase 
space factor 2q/\fs by 2q /So we consider the amplitude 



( 13 ) 
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Figure 4: The upper half-circuit I in the g-plane is equivalent to the hrst full 
circuit I in the s plane, and the lower half-circuit II in q is equivalent to the 
second circuit II in s. 

where = 4 -|- + 47 ^ and 


qn = qo- id- (14) 

This /r(s) is essentially the same as /(s) of (11), but without the ^/s sin¬ 
gularity in p{s). It satishes the unitarity relation (6) with the phase space 
factor 2q/y/s^. In addition to the branch point at s = 4, /r(s) has poles 
at ? = and q = —q^, as shown in hgure 5, both of which have negative 
imaginary parts. It follows that such a resonant amplitude has two poles in 
the second s-sheet. Bearing in mind that the physical region is just above 
the real s > 4 axis in the hrst s-sheet, which is also just above the real q > 0 
axis in q, we see that the pole at ^r in the second s-sheet is near the physical 
region, but the pole at —q^ on the second s-sheet is far from the physical 
region (in the sense of distance travelled in the complex plane). These poles 
in s are at s = 4 -|- 4g|^ (position A in hgure 5) and at s = 4 -|- (position 
B in hgure 5). 

2.1.4 Unitarity and discontinuities 

The functions /(s) and /r(s) (and all amplitudes we deal with) satisfy an 
important condition called Hermitian analyticity: 

ns) = f{s*)- ( 15 ) 
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Figure 5: The resonance poles in the lower half g-plane both correspond to 
poles in sheet II in the s-plane; only the pole at position A is close to the 
physical region. 
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Figure 6: The discontinuity relation (19). 


Consider for example the function /r(s), and take a point just above the cut 
at position s + ie. Here the square root function (on the hrst sheet) takes the 
value +\/s — 4 as e —)■ 0, and 


/r(s + ie) 


7 \/^ 

Qo + 1^ - - 2iq' 


(16) 


On the other hand, /r evaluated at the complex conjugate position s — ie is 


/r(s 


le = 


go + 7^ - + 2i7g ’ 


(17) 


because, as we saw, the square root function takes the value —y /s — 4 at a 
point just under the cut. So clearly 


/r(s - ie) = /r(s + ie) 


(18) 


and (15) is satisfied. 

Applied to the amplitude T{s) of section 2.1.1, the Hermitian analyticity 
condition allows us to rewrite the unitarity condition (5) as 

r+ - T_ = 2ipT+T_ (19) 

where T± = T(s ± ie), and p on the RHS of (19) is understood to be p+. 
The LHS of (19) is the difference between the values of T just above and just 
below the s > 4 cut: it is the discontinuity of T across the cut. Rewriting 
unitarity equations as discontinuity relations will be an essential tool when 
we come to combine unitarity with analyticity by writing dispersion relations 
for our amplitudes. We represent (19) diagrammatically by hgure 6, where 
the lines with dots on are “on-shell” - i.e. they are physical intermediate 
state particles, not Feynman propagators. 
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2.1.5 The i^-matrix 

Dividing both sides of (19) by T+T_ we find 

T-i _ j,-i ^ _2j^ ('20) 

which is another way of writing the unitarity condition. Since, as we have 
seen, 

P+- P- = 2p (21) 

(20) may be satished by simply writing 

T-^ = K-^-ip (22) 

where K~^ has no branch point (is a regular function) at s = 4. Comparing 
(22) with (10) we can identify K~^ with pcotS. And if we choose 

K = gVisn - s) (23) 

we recover the B-W amplitude (11). Equally, we can take K = 7v^/(5'o + 
7 ^ — g^) with p = ‘Iqj and recover /r(s). In general, 

T = K{\- \pK)-^ = (1 - \Kp)-^K. (24) 

T will satisfy (4) if K is real. In the case of a resonance, we may think of 
K as representing a bound state, coupling to the initial and hnal states with 
coupling g, the factor (1 — ipK)“^ then accounting for the state’s decay to 
the open 2-body channel. 

It is important to note that while (21) is certainly true, the function p{s) 
is by no means the only one that has the required discontinuity 2p. In section 
4.1 we will see how to manufacture a function that has this discontinuity but 
does not have the kinematical singularity at s = 0. 


2.2 Several channels and resonances 


Suppose now that we have two resonances, but still only one channel. We 
might think of adding two B-Ws together to form the amplitude 


fis) 


qI ^ gl 

si-s- \pgl S 2 - s- ipgl' 


(25) 
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But you can soon convince yourself that this will not satisfy the unitarity 
constraint (4). A little more work shows that the violations of (4) are of 
order \g^/{si — <S 2 )|, where g is the larger of gi and g 2 . So if the resonances 
are narrow and well separated, adding the two B-Ws will be a reasonable 
approximation. In cases where the resonances have more overlap, we can 
ensure unitarity by putting the two states into K and letting the machinery 
see to unitarity: 

K = ; (26) 

Si — S S2 — S 

inserting this into (24) guarantees a unitary T. 

This formalism really shows its usefulness when more than one channel is 
open. The quantities K and T now become matrices in the space of channels, 
and so does p which is a diagonal matrix of the form 

/ Pi 0 0 \ 

0 P 2 0 (27) 

V 0 0 P3 / 


in a 3-channel case, for example, with pi ~ [s — (mi -|- m 2 )^]^^^ and p 2 ~ 
[s — (m 3 -|-m 4 )^]^/^, etc. Equation (4) still holds, with T now a matrix , as do 
(22) and (24). The matrix K is real, and it can be shown that time reversal 
invariance requires it to be symmetric [3]. 

For example, in the case of a 1-resonance 2-channel problem, we would 
set 

(qj = l, 2 ) (28) 

Sk- s 


where the p* represents the coupling of the resonance to channel i. Then we 
hnd, for example. 


rr _ 9 i92 

-^12 *2 *2 

- 5 - igtpi - ig^p2 


(29) 


which is often called the Flatte form [4]. Notice that each of the thresholds in 
Pi and P 2 generates two Reimann sheets, so such an amplitude takes values 
on four Reimann sheets. 

We can equally easily deal with the case of more than one resonance in 
more than one channel. For two resonances, we simply set 




PiaPaj ^ QibQbj 
Sa S Si) S 


( 30 ) 


and crank the T = (1 — iA"p) handle. 
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Figure 7: Final state interaction amplitude F{s). 



Figure 8: Production process leading to final s-channel. 


Thus far we have chosen K to represent only one or more resonances. 
There is nothing to stop us including a non-resonant “background” term in 
K, which can be any real function of s without the unitarity-induced branch 
points. 

3 Unitarity in Two-Hadron Final State Inter¬ 
actions 

We now consider an amplitude F{s) represented by figure 7, where there 
is for the moment only one final state channel, and where the wiggly line 
could stand for a one-particle state, or for the partial wave projection of 
a two-particle state amplitude, or for the projection of a more complicated 
production amplitude such as the one shown in figure 8, in which we are 
going to parametrise the blob as just some “production vertex”. 

Consider for example a decay of the form K —)■ tt -|- tt, where the charge 
labels are irrelevant for the present purpose. We can picture this as proceed¬ 
ing via an initial weak transition to a two-pion state, followed by recattering 
of the pions in the final state, because mx is above the two-pion continuum 
threshold at 2m.,^. In this case F{s) will be the full decay including the strong 
rescatterings, and T(s) will be the elastic two-body tttt Tin amplitude de- 
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Figure 9: Discontinuity relation for F{s). 


scribing the rescatterings. F{s) will here be evaluated at the discrete point 
s = but in a process such as that in hgure 6, the hnal state variable s 
will run continuously over a phase space interval. 

Let’s suppose that the two final state particles scatter via a strong inter¬ 
action T-matrix T{s). Then just as in (19) and figure 6, F{s) will satisfy a 
discontinuity relation 


F+ — FL = 2iT+pF_ = 2iT_pF+ 
which is equivalent to the unitarity constraint 


(31) 


ImF = TpF* = T*pF. (32) 

Relation (31) is represented by figure 9. Since ImF must be real, (32) shows 
that F must have the phase of T. This is an important result, known as 
Watson’s theorem [5]. 

We might suspect that there is a “F-matrix” type of solution to the 
unitarity constraint on F. Indeed, writing (31) 


(1 — 2iT_p)F+ = F_ 


(33) 


and substituting for T_ from (24) we find 


1 -|- iFp 

which by inspection is satisfied by 


(1 - iA»F+ = 


Fx = 


1 — iKp 


-P 


(34) 


(35) 


where P has no branch point at p = 0. 

How would we describe the production of a single resonance in this for¬ 
malism? We know that for K we would take K = p^/(sr — s). If we were to 
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choose just a constant, say, for P, we would end up with a zero in F at the 
(supposed) peak position sr. Instead, we should take 


fp9 

Sr- s 


(36) 


where fp represents the coupling of the resonance to the initial state, and g 
is as before its coupling to the hnal 2-particle state. A verihcation that (36) 
is the right prescription can be provided via potential theory [6]. Then of 
course we obtain 


F 


fp9 

Sr - S - ipg"^ ■ 


(37) 


The generalisation to the multi-channel case is straightforward: K re¬ 
mains the same matrix as in the elastic multichannel 2 —)■ 2 case, and P 
becomes a column vector with a single channel index, since it describes pro¬ 
duction from a fixed initial state p to a variety of hnal states i. Thus for the 
production of several resonances a = 1, 2 ,... decaying to several hnal states 
i we take 

Pip = Eagia -^- fap (38) 

Sa- S 


and P = (1 — iKp)~^P as a matrix equation. 

All this is quite simple - but, less obviously, we can include background 
terms in both K and P and still be sure that the end result obeys unitarity. 
Suppose for example we take 


P = + P, K = (39) 

Sr — S Sr — S 

in a one-channel problem (with an obvious multi-channel extension). Then 

p =-^ ~ = {fp/g)P^ sin 5 + cos 5. (40) 

Sr- s- ig^p Sr - s - ig^p 


If B is real, the phase of P is still that of the elastic scattering amplitude, 
but there is now a “cos S” piece, as well as a piece proportional to the BW 
amplitude. However there is no requirement that B actually has to be real: 
it might, for instance, be taken to be a Deck amplitude [7], which involves 
a complex Reggeised pion exchange production process. This kind of model 
was used [8] in an early analysis of dihractive production of the ai (one pole, 
one channel), the Pi(1270) and ^2(1400) (two poles, two channels), and 
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some N* states (one pole, two channels). A similar but more elaborate two 
channel analysis of the Oi was done by Basdevant and Berger [9]. 

One can also add a term of the form aBK to P, which would represent 
production by the process B followed by elastic scattering via T - i.e. a term 
aPexpi^sin^ in (40). This was done in the analysis of the ACCMOR data 
by Daum et al. [10]; they concluded that the behaviour of the pn L = 0 

= 1 + amplitude could be explained in detail by such a model, in which 
the Deck amplitude is rescattered through the ai, which may also be directly 
produced. 

A particularly good (and more recent) example of P-matrix methods is 
the = (OO)’''^ state of 7 r’'' 7 r“. Anisovich and Sarantsev [ 11 ] made a 

£t to the available scattering date from nn threshold up to i/i = 1.9 GeV. 
They included 5 channels {i = tttt, KK, dvr, 7777 and 7777 ') and 5 poles in their K- 
matrix. They also added a slowly varying term, and an Adler zero term. The 
FOCUS collaboration [ 12 ] used this iF-matrix and formula (35) to describe 
this TTTT wave in a Dalitz plot analysis of D decays to (only the hnal 

i = TTTT channel is required from F). The P vector contained the same 5 
poles, with 5 new coupling parameters fap, and a new background term (but 
no Adler zero). The collaboration found, in particular, that the low mass 
7 r''' 7 r“ structure of the D~^ Dalitz plot was well reproduced in this iF-matrix 
model, without the need of an ad hoc “a” term. This shows that the same 
7F-matrix description gives a coherent picture of both two-body scattering 
experiments involving light quark constituents and (heavier quark) charm 
meson decays. 

A second interesting application of this formalism is to the Ktiti system 
in B decays to iFi(1270)7r and iF2(1400)7r, where the Ks have = 1+ 
and decay to Ktitt [13]. Following the analysis by the ACCMOR collab¬ 
oration [14], the Kttti system is described by a iF-matrix model compris¬ 
ing six channels (iF*(892)7r, pK, iFg(1430)7r, /o(1370)iF, (iF*(892)7r) 1 = 2 , ^K)-, 
and two resonances. (We note in parenthesis that these channels (and the 
ealier up one) are only “quasi two-body” channels, since one of the two parti¬ 
cles is strongly unstable; we shall discuss this point further in section 8 ). The 
/F-matrix parameters were determined from an ab initio £t to the ACCMOR 
data from the WA3 experiment {K~p —)■ K~'k^'k~p at 63 GeV), taking 

P={1+tK)B + R (41) 

where R is the two-pole direct production term having the form (30), and B 
is a background term (as usual (41) is understood to be a vector in channel 
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Figure 10: Cauchy’s theorem contour. 

space). This will generate an amplitude F via (35). In applying the model to 
B decays, the same iC-matrix was retained (as in the FOCUS work), but the 
background B was set to zero, since no strong diffractive production process 
was now present. 

So far we have been considering parametrisations of the data that re¬ 
spected the constraints of two-body unitarity. We turn now to the inclusion 
of another ingredient - analyticity. 


4 Combining Unitarity and Analyticity 


4.1 Elastic two hadron —> two hadron reactions 

4.1.1 The Chew-Mandelstam phase space factor 


Let’s begin by briefly recalling some simple formulae. Suppose f{z) is ana¬ 
lytic in and on a closed contour C (see hgure 10). Then Cauchy’s theorem 
implies that f{z) can be written as 


fU) 


1- f ibid. 

27ri Jc' z' — z 


(42) 


Now suppose that f{z) has only one branch point at the real value s = 4, 
with a cut attached running along the real axis, s > 4. Then we can freely 
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Figure 11: The distorted contour C, consisting of a circle at infinity and a 
piece from +oo to 4 going below the cut, plus a piece from 4 to +oo going 
above the cut. 


distort C", without running into any singularity of f{z'), into the contour C 
shown in hgure 11. The representation (42) now becomes 


fiz) 


2711 ^4 s ' — Z 


(43) 


assuming that convergence is such that we can throw away the part of C 
at inhnity. The numerator of the integrand in (43) is precisely what we 
have been calling /+ — /_, the discontinuity across the s > 4 cut, which 
is determined by unitarity. So perhaps we can construct a more complete 
parametrisation by combining unitarity with analyticity. 

In fact, equation (20) tells us that the discontinuity of the inverse of the 
elastic amplitude is determined only by the phase space factor. Let us see 
where this leads us. Applying (43) to T“^(s) we obtain 



I{s). 


(44) 
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Unfortunately the integral diverges logarithmically, but if we are content to 
input an arbitrary constant into the calculation in the form of the value of 
T~^ at s = So we can write 


T-\s) 


+ (-^(s) - I{so)) 




So) 


TT 


s' -4 


ds' 


s' (s'- So) (s'-s)’ 


(45) 

(46) 


where the integral now converges. It is convenient to take So = 4. Then we 
hnd 

T“^(s) = constant + L{s) (47) 


where 


L(s) 




~\/ s — 4 + i/s/ 


(48) 


and the logarithm is dehned so that its imaginary part is —tt for s real and 
greater than 4. A careful study of L(s) shows that despite appearances it 
does not have the branch point at s = 0 present in p(s). The imaginary 
part of the logarithm correctly reproduces the unitarity requirement, and 
the integration in I (s) has banished the singularity at s = 0 to an unphysical 
sheet. Functions such as L{s) were introduced by Chew and Mandelstam 
[15]. Note that more generally we would still satisfy unitarity if we replaced 
“constant” in (47) by a regular function r{s) where r has no RH cut and can 
be identihed with K~^. This produces a T{s) whose only branch point is at 
s = 4, but it may of course have resonance poles on the second sheet reached 
through the s > 4 cut. 


4.1.2 Reconstructing the resonance amplitude /r(s) 

A somewhat more complicated exercise in the use of (43) is provided by the 
resonance amplitude /r(<s) of (13). We have 


fn{s + ie) = (7yg^/2go) I ^ ^ - - — I, ga = go - U- (49) 

I g + gR g - gn J 

Along the lower side of the cut, g is replaced by —g so that the discontinuity 
of /r is 


/r+ - /r- = (7y^/2go) 


-1-1—I- 

g + gR g - gR -g + g^ -g - gR 
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where 


(g2 


-qIW - 


- 2igo7 


and 


S(g) = 


272^5; g 




( g 2 - g 2)(^2 

Then according to (43) it should be the case that 


/r(s + ie) = 


TT 


dg' 


n 


£(«') 


g'2 — g2 — ie 


(50) 

(51) 

(52) 


(53) 


The reader may verify by contour integration that the right hand side of (53) 
does indeed reconstruct /r. 

The integration in (53) is understood to be along the top side of the s' > 4 
cut. It will be convenient in section 5.4.2 to consider the integration to be 
running just below the cut instead. In that case, g' will be replaced by —g' 
and we will have the representation 


1 /‘°° 

/r(^+) = - / 

^ J 4,below 


, S(i»') 

s — s' + ie 


(54) 


where 

(S-/2)(S- J*2) 

and 

/2 = 4 + 4g2 . (56) 



The contour for (54) is shown in hgure 12, which also exhibits the poles of 
the discontinuity function S(s'). 


4.1.3 Left hand cut, the N and D functions 


Actually, amplitudes have other singularities in addition to those generated 
by unitarity - in particular, partial wave amplitudes have “left hand” sin¬ 
gularities associated with exchange processes, which typically produce cuts 
along the real s axis for s < sl, where sr < 4. We can then write 


T(s) = 



D{s) 


(57) 
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A2 = 4 


A2 


j*2 



Figure 12: The contour in the A^-plane for the integration in (54). 


where N{s) has only the LH cut and D{s) has only the RH cut. In that case, 
D+-D_ = N{T-^ - T:^) = -2ipN (58) 

and we can therefore write (compare (44)), 


TT 44 s' - S 


(59) 


assuming that N{s) is such that the integral converges. 

Another representation for D is sometimes useful. From (58) we obtain 


D+ — D- = —2i H+e'^ sin d (60) 

using (57), so that 

n+ = = D_S-^ (61) 

where S is the 1 x 1 S' - matrix: S' = 1 + 2iT p. Taking the logarithm of the 
first equation in (61), we see that the discontinuity of \nD is —2iS, and so 
we can write 

1 r°° S(s') 

D{s)=exp{— ^ds'}, (62) 

TT 44 s' - S 

always assuming convergence. 

We shall not pursue the calculation of 2 —)■ 2 amplitudes any further here. 
Instead our aim will be to see how unitarity and analyticity can provide useful 
formulae for the analysis of hadronic hnal state interactions, going beyond 
the simple 4F-matrix methods so far discussed. 
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4.2 Two hadron final state interactions 

Let’s return to the one-channel hnal state interaction (f.s.i.) discontinuity 
relation, written out again, 

F+-F_= 2iT_pF+ = 2iT+pF_, (63) 

which we previously arranged to satisfy in the K - matrix / P - vector 
formalism. This time we’re gong to include analyticity, as in our discussion 
of L{s). 

First note that we can write (63) as 

F+ = (1 + 2iT+p)F_ = SF_ (64) 

But we also know from (61) that S = D_/D^. It follows that 

F+D+ -F_D_= 0, (65) 

so that the function FD has no branch point at s = 4. Hence the unitarity 
constraint is satisfied by 


F{s) = C{s)/D{s) 


( 66 ) 


where C'(s) is any function regular at s = 4, for example a polynomial. 

Now suppose that F{s) has a “background” term B{s) which we want to 
include - for instance, a Deck-type production process; B{s) is assumed to 
have only a LH cut. We would like to take account of both B{s) and the 
unitarity constraint, in a way consistent with analyticity. In this case we can 
satisfy our unitarity and analyticity constraints by writing 


F{s) = B{s) + 


r(s')p(s')F(8') 


TT 


ds'. 


(67) 


which is an integral equation for F{s). Remarkably, there is an exact solution 
of this equation, due to Omnes [16] and Muskhelishvili [17]. 

Consider the discontinuity of the quantity D{F — B) across the elastic 
s > 4 cut: 


disc(DF - DB) = D+F+ - D_F_ - (D+ - D_)B 

= -{D+ - D_)B = 2ipNB. (68) 
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We can therefore write 


D+{F+-B) = 


TT 


p{s')N{s')B{s') 
s' — s — ie 


ds' 


or 


p'N'B' 


F{s) = B 


ttD 


ds'. 


(69) 

(70) 


This is the famous 0-M solution to our f.s.i. problem. 

There are three points to note immediately: 

1) If i? is a constant (and the integral in (70) converges), then F becomes 
simply B/D. 

2) We can always add to (70) any solution to the “homogeneous” version 
of (67) - that is, the equation with B{s) = 0. We already know that such a 
solution has the form F{s) = C{s)/D{s) with C(s) regular at s = 4. So we 
may write the general solution to our problem as 


F{s) = B{s) + 


TiD{s) Ji : 
3) By making use of the identity 

1 P.V. 




Ie 


D{s) 


(71) 


s' — s — le 


s' — s 


+ i7r(5(s' — s) 


(72) 


where “P.V.” stands for “principal part”, and taking N to be approximately 
constant, we can write (71) as 


F BF^ cos 5 + 



- p y roo 

itH-/ —;- 

L JA s' - S \ 


(73) 


which is the more sophisticated version of our iP-matrix formula that analyt- 
icity has bought us. The task of htting (73) to data is simplihed by the fact 
that the principal value integral is independent of the 2-body scattering pa¬ 
rameters. Expression (58) was used by Bowler et al. [8] to model diffractive 
ai production, using a Deck amplitude for B{s). 

As with the iP-matrix approach, the foregoing can be extended to the case 
of several coupled 2-body channels, so that T, N, D and p become matrices 
in channel space, while F and B are vectors. We shall not discuss 2-body 
hnal states further here, but turn our attention now to our main topic, the 
problem of f.s.i. among three hadrons. 
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Pi 


P2 


P3 


Figure 13: Three-hadron f.s.i. amplitude. 


5 Final State Interactions Among Three Hadrons 

5.1 Kinematics, and the isobar model 

The amplitude we are now concerned with can be represented by hgure 10, 
in which a state of dehnite decays to three hadrons. We hrst have to 
understand the kinematics of such decays. To simplify matters, we shall 
suppose that the three final state particles are spinless and of equal unit 
mass, while the initial state has = O’*" and invariant mass m (i.e. this is 
the energy in the 3-body c.m. frame). 

Let the 4-momenta of the final state particles be pi,P 2 and pa, and let P 
be that of the initial state, so that P = Pi + P 2 + Ps, with P^ = rrp. We 
introduce invariant variables s, t, u by 



(74) 


(75) 



2 p(s, rrP)q[s)Xs 


(76) 


where Xg is the cosine of the angle between and Pg in this system, and 
where 


p{s,rrP) = {[s — (m — l)^][s — {m + l)^]}^/^/2v/s = k{s,frp)/2y/s (77) 


,(s) = (s-4)'/V2. 


(78) 
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So q{s) is the magnitude of the momentum of particle 2 or 3 in the 2-3 
c.m.s., and p{s,m?) is the magnitude of the momentum of particle 1 in this 
system. The physical region for the decay process m—)-l-|-2-|-3is then 
4 < s < (m — 1)^ and < 1; the second condition can be written, after 
some algebra, as 

stu — {w? — 1)^ > 0 (79) 

or, using (75) as 

r(s, t, m^) = st{3 + 771^ — s — t) — {m^ — 1)^ > 0. (80) 


r is the Kibble [18] cubic, drawn in hgure 14 for the case m > 3. The physical 
region for the decay (the Dalitz region) is inside the closed loop labelled V. 
The regions labelled I, II and II are physical regions for the ‘crossed’ reactions 
m + 1 ^ 2 + 3 (I), m + 2 ^ 1 + 3 (II), m + 3 ^ 1 + 2 (III). 

A hgure showing how experimental events populate the Dalitz region is 
called a Dalitz plot - hrst invented by Dalitz [19] in connection with his 
famous analysis of 77 —)■ tttttt (which led to the discovery of parity viola¬ 
tion). If the matrix element for the decay is a constant, then the events will 
be uniformly distributed on the plot. This follows from the fact that the 
dependence on s and t in the three-particle phase space is proportional to 

dsdt 
w? 

If, on the other hand, a 2-particle resonance can be formed in any of the 
three pairs of particles, then there will be strong concentrations of events 
along bands centred on the square of the resonance mass(es). 

It is an empirical fact, of course, that very many three-hadron systems 
are such that their two-hadron subsystems do indeed form resonances (e.g. 

, etc.) This fact is the basis for the isobar model [20], 
which expresses the m —)■ 1 -|- 2 -|- 3 decay amplitude as a coherent superposi¬ 
tion of “resonance -|- spectator particle” states, as in figure 15. In practical 
applications, one has to deal seriously with all the complications of spin, an¬ 
gular momentum, recoupling coefficients etc. Eventually we shall come back 
to these complications, briefly, but for the most part we shall for pedagogi¬ 
cal purposes concentrate on the simplest model, in which all three identical 
particles are spinless, each pair forms a resonance in the / = 0 state, and the 
overall = O’*". For this toy isobar model, then, referring to hgure 15 we 
write 



F{s, t, u, m^) = C{m?‘)M{s) + C{m?‘)M{t) + C{m?‘)M{u) (82) 
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t = {m+lf 



Figure 14: The Kibble cubic, with the region for the decay m —)• 1 + 2 + 3 
labelled V (the Dalitz plot region). 




Figure 15: Isobar model expansion of F as a coherent superposition of reso¬ 
nance + particle amplitudes. 
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Figure 16: Subenergy discontinuity relation. 


where C{m^) is a common “production vertex” represented by the solid blob, 
and depends on m? but not on s, while M(s) is the elastic 2 —)■ 2 amplitude 
in the 2+3 channel having only the s > 4 cut (and similarly for M{t) and 
M{u)). The factorisation of each term in ((82)) into a product of a function 
of times a function of s is fundamental to the isobar model - and is also, 
as we shall soon see, inconsistent with unitarity. 

5.2 The isobar model violates subenergy unitarity 

The constraint which determines the correct phase in a subenergy channel 
is the subenergy discontinuity (unitarity) relation, shown diagrammatically 
in hgure 16 for the s channel. This is very similar to figure 9, but with one 
crucial difference: the amplitude F{s,t,u,m‘^) now contains (unlike F{s) of 
figure 9) an angle dependence via (76), which must be integrated over in the 
two-particle phase space integral: 


s > 4 : F{s+, t, u, m^) — F{s_,t, u, m^) = 


2ip(s) 


da;sF(s+, t(s+, x*), m(s+, x*), m^)M(s_). 


(83) 


'-1 
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Equation (83) is the required generalisation of (31). Recall that our normal¬ 
ization for p{s) is 

Pis) = ■ (84) 

It is quite simple to see that (82) cannot satisfy (83). By construction, 
M{t) and M{u) have only the cuts t > 4, u > 4, and no discontinuity in s. 
So the LHS of (83) is just 

= C{m^)2ip{s)M+{s)M_{s). (85) 

This is the same as the RHS of (83) only if the M{t) and M{u) terms in F 
are absent. 

The expansion (82) therefore has to be modihed in order for F to satisfy 
(83). An economical way to do this is to replace M{s) in (82) by M{s)(f){s, m?) 
(and similarly for M{t) and M{u)), where m?) is going to be determined 
from (83) and analyticity. We have anticipated, and will soon conhrm, that 
the correction function 0 will depend on both s and w?, and will thus spoil 
the factorisation in (82) which was alluded to earlier. 

So we now take 

F{s,t,u,m^) = C{m^)[M{s)(j){s,m^) + + M{u)(l){u,m^)]. (86) 

Inserting (86) into (83) the LHS is (in shortened notation) 

C'(M+0+- M_0_) = C'((M+-M_)0+ + M_(0+-0_)) 

= C'(2ipM+M_0+ + M_(0+-0_)) (87) 

while the RHS is 

1 

2\pCM+M_ (j)++ 2\p2C M_ - / M{t)(j){t,m‘^)dxs (88) 

2 

since the two contributions from the t and u terms are equal. It follows that 

s > 4 : 0(s+,m^) — (t){s-^rn?‘) = 2ip(s+) j M(t)0(f,m^)da;<i. (89) 

This important equation tells us that the function 0(s, m^) has a discontinu¬ 
ity across the s > 4 cut, determined to be (89) by unitarity in the s-channel. 
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It is evident from (89) that 0 must depend on as well as on s, via the 
dependence of t (c.f. (76)). Equation (89) also implies that 0 will develop 
an imaginary part for s > 4, which means that the phases carried by the 
terms in (86) are no longer those carried by the two-body amplitudes M 
in (82). The importance of examining the constraint of subenergy unitarity 
was especially emphasized by Aaron and Amado [21]. 


5.3 Implementing subenergy unitarity and analyticity 

5.3.1 Failure of a “A'-matrix” approach 


We might be tempted to implement the subenergy unitarity constraint by a 
A'-matrix type of procedure, exploiting the fact (as before) that p+—p_ = 2p. 
We would then write 


? 


= (regular function) + ip M{t)(j){t,m^)dxs 


(90) 




Unfortunately, though simple, this prescription is incorrect, in the sense that 
it leads to singularities in 0 which perturbation theory teaches us should not 
be impacting the physical region. To see this, note first that the uncorrected 
isobar model corresponds to taking 0=1. So we will make (90) more precise 
by writing 

? 


0(s, m^) = 1-I-ip / M{t)(t){t^m^)dxs 


(91) 




and imagine solving this integral equation iteratively. The first correction 
will be 


J(s,m^)= ip / M{t{s, Xs))dXs. 


(92) 




Suppose that M has one resonance, which for the moment we may parametrise 


as 


M{t{s,Xs)) = 


47\/57 


t2 2 

/ =mR 


iT. 


(93) 


P -t{s,Xs)' 

The denominator is a linear function of Xg, and the integral is easily done 
yielding the result J{s,riP) = 2ip{s)B{s,m?) where 


B{s,m‘^) = 


47 ys; 


In 


m 


R 


t{s,Xs = -fl) - iT 


4p(s,m^)g(s) — t(s, = —1) — iT 


(94) 


^This form ignores the normal threshold branch point at t = 4 present in /(t) given by 
(11) or in /R(t) of (13). We will treat this properly in section 5.4. 
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Figure 17: First rescattering correction, or triangle graph. 


The logarithm has singularities in s when 

t{s, Xs = ±1) = — iF. (95) 

The two curves t{s,Xs = +1) and t{s,Xs = —1) together form the boundary 
of the s — t Dalitz plot as s varies. The singularity at t{s, Xg = +1) = —iF 

occurs at that value of s, say Sh, at which the boundary arc t{s, Xg = +1) hits 
the resonance band centred at f similarly for the other singularity, at 

Though “only” logarithms, these singularities cause quite noticeable phase 
and modulus variation - but are they to be believed? The answer is no: these 
singularities are largely spurious [22]. They are in fact the positions of singu¬ 
larities (in s) of the triangle graph shown in hgure 17. Such diagrams have 
been thoroughly discussed (and we shall soon meet them). One singularity, 
Sb, is near threshold and can be near the physical sheet of cj), but its proximity 
to threshold limits its effect. The other, at is far from the physical region. 

5.3.2 Combining unitarity and analyticity 

Why are singularities openly present in B{s,m‘^) somehow masked in the 
triangle graph? Just as in the case of the s = 0 singularity of p(s), we need 
to include analyticity, in addition to unitarity, in order to get a physically 
correct amplitude. In other words, we have to insert the discontinuity relation 
(89) into a dispersion relation. Singularities present in the integrand can get 
moved away from the physical region after integration. 
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This leads immediately to the equation 

4>{s,m^) = 1 + - [ _ {1 f dxsM{t{s\xs))(p{t{s',Xs),m^)}ds' 

(96) 

where we are assuming that the integral over s' will converge, and are also 
taking the inhomogeneous term to be unity, corresponding to the unmodihed 
isobar model, as in (91). It will sometimes be convenient to define <h(s, m^) = 
M(s)0(s) and rewrite (96) as 

<h(s,m^) = M(s) + 2M(s)- f ^ P^^^ _ {1 f (|)(t(s',a;5))da;4ds'. (97) 

TT J 4 S — S — le 2, J_i 

Equations (96) and (97) are integral equations embodying the basic con¬ 
straints of two-body unitarity and analyticity. They are therefore a kind of 
minimal theory of corrections to the isobar model, which would correspond 
to just the first term in (96) and (97). Equations of this type were first 
proposed by Khuri and Treiman [23], from a rather different standpoint. 


5.4 The first rescattering correction: the triangle graph 

5.4.1 The rescattering amplitude T{s,m?) 

It is possible to proceed directly on the basis of (97), solving it iteratively, for 
example. The first iteration is again just the usual two-body amplitude used 
in the isobar model, and the first correction to this adds to it the amplitude 


2M(s) 

TT 




M(t(s',a;^))da;J. 


(98) 


A particularly interesting case is that in which the amplitude M{t) is reso¬ 
nant. It is important to make sure that we are getting the sheet structure of 
M correct, so we will set M{t) = /R(t) where /R(t) has the representation 

(54): 


M(i) = 


dA‘ 


E(A2 


^ J 4,below 


t -I- i 


le 


(99) 


and 


E(A^) 


167 ^ 1 / 57 a/s — 4 

(A2-/2)(A2-/*2) 


( 100 ) 
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Figure 18: The triangle graph amplitude of (105). 


with 

J2 =4 + 4(g2_^2_2i^^^)_ 

The hrst iteration of (97) then gives 

m?) = M{s)[l + T(s, m^)] 

where 


T(s, m^) = — ds'— 
vr Ja 


, Pis') ir\.ir 


s' — s — ie 2 


dV 


'-1 


TT 


'4,below 


f — + i 


le 


( 101 ) 

( 102 ) 

(103) 


is the first rescattering correction. We have dropped the factor of 2 since we 
want the contribution from just one rescattering channel, and we have also 
suppressed the (important) arguments s',Xs of t. We may rewrite (103) as 


1 r°° 

T(s,m2) = — / S(A2)/(s,A2,m2) (104) 

^ J4,below 


where 


/(5, A 


2 2 \ 

,m ) 


ifd.- i 

TT J 4 s' - s - le 2 



dxs 


1 

A2 — t — ie 


(105) 


is the triangle graph of hgure 18 with two internal particles of unit mass and 
one of squared mass A^. 

We can understand this graphical interpretation by considering how we 
would calculate this diagram by writing a dispersion relation in the variable 
s. Looking along the s-channel, we see a normal threshold at s = 4, with a 
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Figure 19: Reconstructing figure 18 by dispersing in s. 



Figure 20: The first rescattering correction T{s,m?). 

discontinuity given by cutting the graph across the two internal lines of unit 
mass, which are put on mass-shell. This discontinuity is proportional to the 
product of the two-body phase space factor p(s) and the s-wave projection 
(in the s- channel c.m.s.) of the t-channel exchange diagram, as shown in 
hgure 19. These are exactly the ingredients of (105). 

The amplitude T(s,m^) of (104) is therefore an integral over the variable 
internal squared mass of the triangle graph, weighted by the spectral 
function Fl(A^). We may represent T(s,m^) by hgure 20. 
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5.4.2 Singularities of T{s,m?) near the physical region 

An important question is whether T{s,m^) as given by (104) has any sin¬ 
gularities in s or which are near the physical region in those variables, 
since they would be likely to cause signihcant variation in the magnitude 
and phase of T(s, w?). The singularities of amplitudes such as T(s, w?) were 
studied in [24], whose analysis we now briefly describe. . 

It is clear first of all that T has the normal threshold branch point at s = 4, 
since this is present in /(s, As usual, this results in a two-sheeted 

structure for T. The physical amplitude T{s,im?) is obtained from (104) by 
integrating, with s approaching the real axis from above, the physical sheet 
amplitude of /(s, along a contour taken just below the real axis, 

as shown in figure 12. Also in this figure we have indicated the positions of 
the poles in S(A^) at A^ = P and A^ = 

There is another way in which a singularity of T(s, m?) can be generated. 
It can be shown [25] that /(s, A^, w?) has two singularities at A^ = A^(s, rn?) 
which move around in the A^-plane as s (or rp) move. It may happen that, 
as s moves in the complex s-plane, one of these singularities - say A^(s, m?) - 
approaches the A^ contour from above, so that the contour has to be deformed 
away from the advancing singularity in order to have a smooth continuation, 
as shown in figure 21. If it should happen that the advancing A^(s,m^) 
actually pins the contour against the pole of S(A^) at A^ = P, so that the 
contour cannot be deformed away, then for that value of s and rp there will 
be a singularity of T{s,m?). This is called a “pinch” singularity, for obvious 
reasons. 

Careful analysis [25] [24] shows that it is possible for such a pinch singu¬ 
larity of T{s,w?) to occur at a point s = Sbijm?), near the physical region 
in s. The singularity, which is logarithmic, is present on the second s-sheet 
of T{s,rP), reached as usual by crossing the real axis from just above the 
s > 4 cut. The imaginary part of Sb is related to that of P, and for a narrow 
resonance Sb will be close to the real axis and therefore near the physical 
region for T{s,rP). 

This situation only arises for a particular range of values of s and P, for 
fixed rp. This range is easily visualized on a Dalitz plot for the variables 
A^ and s. Referring to figure 22, the pole at A^ = P is represented as a 
resonance band at fixed A^. This band intersects the boundary of the plot 
at two points Sb and Sa (assuming the imaginary part of P is small). Both 
of these points are potentially singularities of T(s, rp) - indeed they are just 
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Figure 21: Pinch singularity in T{s,m?). 

the same points as those we encountered using the (incorrect) “iP-matrix” 
unitarisation scheme. Here, analysis shows [24] that Sa is never near the 
physical region and Sh is only near it if P lies in the range 

— 1) < P < {m — 1)^, (106) 

neglecting the complex part of P. The corresponding Sh lies in the range 

4<Sb<(m + l). (107) 

The ranges (106) and (107) are of course given for our current (unit) mass 
values. In general, the resonance band must intersect the Dalitz plot bound¬ 
ary on the upper left hand arc, and the singularity is read off on the s-axis 
from this intersection. 

The upper limit of the range (106) corresponds to the np value nP = 
(/ -|- 1)^, which is just the “normal threshold” for making the “quasi two- 
particle” state consisting of one particle of complex mass I and another 
particle of unit mass. We shall discuss such particle -|- resonance states in 
section 8. It is clear that no nearby singularity associated with the resonance 
can occur unless rp is at least greater than the threshold value (/ -|- 1)^. For 
this rP, Sb{rP) is close to the point m -|- 1, with a small negative imaginary 
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Figure 22: The situation in which the triangle singularity Sb is near the 
physical region. 
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Figure 23: The motion of Sbirn^) as increases from (J +1) to greater than 
2/2 + 1 . 

part related to that of P. As increases, the Dalitz plot grows, and the 
intersection point Sbim^) moves towards the point s = 4, reaching it when 

= 2/2 + 1. Thereafter Sb{m?‘) moves around s = 4 into the upper half 
plane, still on the second s-sheet, but progressively further from the physical 
region. This motion of Sb{m?) is shown in hgure 23. In general, Sb will he 
close to threshold. 

The amplitude T{s,m^) could be evaluated directly from the represen¬ 
tation (104). However, it would be somewhat more intuitive if we could 
somehow extract from the A2 integral the contribution associated with the 
pole at A2 = /2. We would then, up to a constant factor, be dealing with 
/(s,/2,m2), which is the triangle graph with an internal particle of squared 
mass /2, shown in hgure 17. This can easily be done. Recall that the sin¬ 
gularity at Sb arises from a pinch of the A2 contour in hgure 21 between the 
singularity at A2 = A^(s,m2) of and the pole at A2 = P. If 

the A2 contour passed below the pole, A^(s,m2) and the pole would be on 
the same side of the contour, and no pinch would occur. Let us denote the 
amplitude dehned along such a contour by T{s,m?). Then T is free of the 
singularity at Sb- Referring to hgure 24, we see that this second contour (for 
T) is equivalent to the hrst contour (for T) together with a circuit around 
the point P = P. Hence 

T(s, nP) = f{s, rp) -|- —2ttiR /(s, P, rp) (108) 

TT 

where R is the residue of H(A2) at the pole P = P. The function /(s, P, rP) 
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Figure 24: The contours for (108). 


does contain the nearby singularity Sb, while the function T{s,m?) does not. 
We can therefore calculate the effect of the singularity Sb by evaluating the 
quantity 2ii?/(s, A^, m^). Since R = we obtain hnally for the 

singular part of the rescattering amplitude 


^sing ( ^^ 




m 


(109) 


Noting now that — 7 ^ — 2igo7)^'^^ ~ —qo, we see that we have 

arrived at the reassuring result that is, to a good approximation, just 
the triangle graph of figure 17, using the “naive” amplitude 47 y^/(/^ — t) 
for the f-channel resonance (i.e. ignoring the branch point at f = 4): 


T- Um^) « 1 /“ds'—dfl_i /‘TVS 

^ /•OO 


dx 


71 




s' — s — le 


( 110 ) 


where B{s,m?) is given in (94). 


5.4.3 Physical picture of the nearby rescattering singularity 

In the limit where the imaginary part of P goes to zero, so the resonance has 
zero width, the singularity Sb{Tn?) approaches the real axis. This is an infinity 
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in f{s,P,m?) in the physical region (though not of T, in which / is multi¬ 
plied by the width parameter 7). How can such a physical region singularity 
of a Feynman graph occur? The general answer was provided in an elegant 
paper by Coleman and Norton [26]. They showed that a Feynman ampli¬ 
tude has singularities in the physical region if and only if the corresponding 
Feynman diagram can be interpreted as a picture of a four-momentum con¬ 
serving process occurring in space-time, with all internal particles on-shell, 
and moving forward in time. The particular case of this result for the triangle 
diagram was given by Bronzan [27]. 

Following [27] for our simple case of three identical spinless particles of 
unit mass, and a resonance of squared mass /^, consider the rescattering 
graph of figure 17, in which the resonance is in the (13) or t-channel, and the 
hnal rescattering is in the (23) or s-channel. For this to be a real physical 
process, we certainly need tv? > (J -|- 1)^, or /^ < (m — 1)^ as in (106). The 
lower inequality in (106) arises from the “catch-up” condition; namely, in the 
rest frame of P, the decay particle 3 must be moving in the same direction as 
the “fleeing” particle 2, and the speed of particle 3 must be greater than or 
equal to the speed of particle 2. The kinematics is similar to that in section 
5.1, except that now we work in the (13) c.m.s. rather that the (23) c.m.s., 
and we set t = P. So we write (c.f. (76)-(78)) 

s = {?) + riP — t)/2 — 2p(t, rp)q{t)xt (HI) 


where 

p{t,nP) = [t^ — 2t{nP -I- 1) -I- {nP — l)^]^'^^/2\/t (112) 

and 

g(t) = (t-4)1/72, (113) 

and Xt is the cosine of the angle between 1 and 3 in this t-channel c.m.s. 
The speed of particle 3 is then [{P — 4)//^]i/^. The energy of particle 2 is 
{w? — P — l)/2'/P, and the magnitude of its momentum is p{P,rP). We 
therefore require, for a physical rescattering, Xt = +'\- and 


\iP - 4)//7/2 > - 2 / 7 ^^ + 1) + {rP - 1)T^^ 

“ w? — P — 1 


(114) 


which reduces to 

P > {tP - l)/2 (115) 
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Figure 25: The region AB of the Dalitz plot boundary where the rescattering 
can occur as a physical process. 

as in the lower inequality of (106). The arc Xt = +1 is shown in hgure 
25, from which, together with (115) we see that the catch-up conditions are 
precisely that the P band intersects the Dalitz boundary on the upper left- 
hand arc AB. This condition is completely general, for arbitrary mass values 
in the triangle graph. 

5.4.4 Some examples 

Under what circumstances might a nearby singularity si, be potentially ob¬ 
servable? We’ll return to this question in the next section, but hrst we discuss 
some possible examples. 

In the case of three identical hnal state particles, there will be a resonance 
at s = P in M{s), and this amplitude multiplies f{s,P,m‘^) in (58). This 
situation is shown in hgure 26. It is clear that Sb will he far from the region 
where M{s) is large, and the net effect will be only a small modihcation of 
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Figure 26: Identical resonances in the s- and t-channels. 


the tail of the resonance in M{s). 

We need to consider, instead, a case where the s and t channels, say, 
contain different interactions. An example of this type of triangle graph was 
calculated in [28]; see also [29]. The reaction considered was nN —)■ miN, 
and was the triangle shown in hgure 27, where the intermediate 

resonance was the A(1232). All particles were treated as spinless, interacting 
in s-waves only. The triangle was calculated using a dispersion relation in s, 
just as in figure 19 for /(s, see figure 28. Both singularities Sa and 


TT 



Figure 27: A triangle calculated in [28] . 
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Figure 28: Calculating figure 27 by dispersing in s. 

Sh showed up clearly in the integrand of /(s, but only Sb produced 

any effect in |/(s, /^, m?) p, and then only when P was in the expected range 

-(m^ + — 2) < P < {m — 1)^ (US) 

which is the analogue of (106) for this case (the pion mass is unity, and that 
of the nucleon is M). The sharpness of the Sb effect depends sensitively on 
the width of the resonance I. For a realistic A width, no separate peak near 
s = Sb was seen in \f{s,P,m‘^)\^, only a rise near threshold. For a width of 
order one tenth of the true width, a peak in the intensity near Sb was present. 
In practice, since the width is also an overall factor in Tsing(s,m^) of (58), 
there will be a trade-off between the closeness of Sb to the physical region 
and the magnitude of the effect. 

In the 1960s and 1970s considerable effort went into trying to hnd a reac¬ 
tion in which a triangle singularity might be detectable. To my knowledge, 
no such effect was ever conclusively demonstrated in those days. More re¬ 
cently, however, the idea has been revived in various contexts. For example, 
Szczepaniak [30] considers diagrams of the type shown in hgure 29, where 
his notation for the masses is used. He calculates two cases. In the first, 
“M” is the Y(4260) state, “m” is the average of the D and D* masses, “/i” 
is a (massless) pion, and “M'” is the J/V’- For a t-channel resonance Dq at 
mass 2.4 GeV, he finds an enhancement near the s-threshold due to the Sb 
singularity, close to the Zc(3900) seen in the (J/'^tt) hnal state. All spins 
were neglected, and interactions were in s-waves. In a second case, Szczepa¬ 
niak takes “M” to be the T(5s, 1086), “m” to be the average of the B and 
B* masses, and “M'” to be the T(ls). A t-channel B* resonance at 5.698 
GeV produces an enhancement near threshold in the T(ls) 7 r channel, in the 
region of the observed Zc(10610) peak. 


43 







Figure 29: The triangle considered in [30]. 



Figure 30: Case of two different resonances in the s and t channels. 

5.4.5 Enhancements in the three-body (m^) channel 

So far we have concentrated on the possibility of a signihcant effect in the 
subenergy variable s, as varies. In cases where the hnal state allows dif¬ 
ferent resonances in the s- and t-channels, rescatterings of the type shown in 
hgure 30 will occur. Here and are the squared masses of the resonances 
in the two channels. In such a case will effectively be evaluated 

at s = If, and the rescattering amplitude f , m"^) will exhibit a singu¬ 
larity in at = ml say, near to the physical region in if the and 
/| bands cross on the “magic” upper left hand arc of the Dalitz plot (which 
only happens for one value of m?)\ see hgure 31. 

One such process, considered in [28], is shown in hgure 32. In general, 
the m^-enhancement will occur near the threshold, in this case at m? = 
(ruN + fyipY- There are several baryon resonances in this energy region, but 
their dynamical origin is diherent from the triangle singularity. 

Recently it has been suggested [31] that the ai(1420) may be identihed 
with such a triangle enhancement. Here the diagram is shown in hgure 33. 
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Figure 31: The critical value = ml such that the resonances in the s and 
t channels cross on the rescattering arc AB. 


/ TT 





Figure 32: The m^-enhancement triangle considered in [28]. 
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Figure 33: The triangle graph identihed with the ai(1420) in [32], 


Spin and kinematic factors are included, and a peak in |/p, with a sharp 
phase motion relative to a reference wave, is found. The calculated effect 
(roughly 1% in peak intensity) is consistent with the data [32], 


5.4.6 The observability of triangle singularities, and Watson’s 
theorem 


Consider a simple model in which there is a t-channel resonance in the ampli¬ 
tude Mis{t), but not in the s-channel amplitude M 23 (s), and no interaction 
in the u-channeh Then including the first rescattering correction with a 
logarithmic singularity Sf,, we have 


= C'(m^){M 23 (s)[l T(s,m^)] Mi 3 (f)} 


where 


and 


T(s, m^) 





s' — s — le 


5i3(s',m^) = - / da:^Mi3(t) 


'-1 


1 r j 

- / dXo^^-7-r. 

2 J_i /? - t(s', X,) 


(117) 

(118) 


(119) 


The discussion of Watson’s theorem refers to a particular partial wave 
in a two-body channel. So consider the I = 0 partial wave projection of 
F{s,t,m‘^) in the 2-3 cms. This is 

Fo(s,m^) = 2 / F{s,t,m‘^)dxs = C{rri^){M23{s)[l+T{s,m‘^)]+Bi3{s,m‘^)}. 

( 120 ) 
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We now observe that the same function i?i 3 (s,m^) appears in the projection 
(120) and in the dispersion integral for T{s, w?). But whereas - 813 ( 5 , m?) has 
both singularities Sa and near the physical region, as we saw in section 
5.3.1, for T( 5 , m?) only the singularity at Sh may be near the physical region, 
and then only for a range of If. Focusing then on the case in which Sb is 
close to the physical region, we may ask: what is the net effect of having the 
Sb singularity present in both T^s^rri^) and - 813 ( 5 , 771 ^)? 

This was the question raised and answered by Schmid [33], and further 
discussed in [34]. Using the identity (72), we can always write T(s,m^) as 


T{s,m^) = ip{s)Bi 3 {s,m^) + 


P.V. 


71 


ds 




s' — s 


( 121 ) 


The surprising fact is that it can be shown [33] [34] that near the point 
s = Sb, when Sb is near the physical region, the Principal Value integral in 
(121) contributes equally with the (5-function, so that 

T(s ^ Sb,m^) ^ 2ip(s)Bi3(s,m^). (122) 


Then FQ{s,m?) of (120) becomes 

Fo{s,m^) K. C{m^){M2z{s) + Bi-i{s,m^)[l + 2ip{s)M2z{s)]} 

= C'(m2)[M23(5) +- 813 ( 5 , m2)e2‘^^3(")]. (123) 


Thus the net effect of the nearby singularity Sb in the rescattering correction 
to the projected amplitude Fq is simply to modify the phase of the projec¬ 
tion, - 813 ( 5 , 771 ^), of the t-channel resonance. This was Schmid’s result [33], 
conhrmed in [34]. 

Put differently, the intensity without the rescattering would be propor¬ 
tional to |M23(5) -|- - 813 ( 5 , 777 ^) 1 ^, and with the rescattering to |M23(5) -|- 
g2i523b)^^g('_g^ 777 ^) p. While these may differ in magnitude, the presence of the 
rescattering singularity in T{s, m?) cannot be distinguished from its presence 
in - 813 ( 5 , 777 ^). It would seem that the only surviving observable effect of the 
triangle singularity is the modihcation of the interference between M 23 {s) 
and - 813 ( 5 , 777 ^). Though a subtle effect, it may be relevant to experiments 
seeking to extract phase information from Dalitz plot interferences. 

In concluding this section, we return to Watson’s theorem. It is clear 
from (58) that the phase of ^ 0 ( 5 , 777 ^) is certainly not ( 523 ( 5 ) for two reasons: 
hrst, the projection - 813 ( 5 , 777 ^) is complex, and second so is the rescattering 
term T{s, rn?). 
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5.5 The single variable representation for $ (1) 

Although, as we said, one could proceed on the basis of equations like (97) 
as it stands, this involves a double integral on the RHS. It seems desirable 
to convert (97) into a single variable integral equation, if possible. This will 
be more convenient for numerical work, and it also turns out to be much 
better suited for discussing general properties of the model - in particular 
the perhaps surprising fact that it can satisfy three-body unitarity as well. So 
we now turn to the single variable representation [35] for $ (or 0). 

The single variable representation (SVR) was first obtained in [35], and 
we shall outline that derivation here. In section 5.5 we shall discuss an 
alternative, more general, derivation given by Pasquier and Pasquier [36]. 

In the present approach, the key step (due to Anisovich [37]) exploits the 
fact that <h(f, m^) is analytic in the t-plane cut along the real axis f > 4, so 
that we can write (always assuming convergence) 


<h(t, m^) 


1 

27ri 



<h(A^, m?) 
A2-t 


dA^ 


(124) 


where the contour C loops around the cut in a clockwise sense (c.f. figure 
11). Then the integral term in (97) becomes 


2 


M(g) 

27r^i 


r p{s') 1 

^4 s' - s - ie 2 






(125) 


which seems to have made matters worse. But we are going to invert the 
orders of integration in (125), after which things will look better. 

To do that, we need to be careful about the way the various singulari¬ 
ties of the integrand are situated, with respect to the integration contours. 
A useful trick is to use a form of the third ingredient of A-matrix theory, 
namely crossing symmetry. In the present case, this will assert that our de¬ 
cay amplitude for m —)■ 1 -|- 2 -|- 3, with m > 3, is the analytic continuation 
in m of the 2 —)■ 2 amplitude for m -|- I —)■ 2 -|- 3 which is shown in figure 
16. In practice, this means starting at a value m < 3, where the decay is 
not possible, doing the contour shuffling, and continuing the result to a value 
m> 3. After this manoevre, (125) becomes 

[ $(A^m^)/(s, A^m^)dA^ (126) 

J c 
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Figure 34: The crossed process m + 1 —)■ 2 + 3. 


which indeed has <F under a single integral, multiplied by the function 


/(s, A^m^) = — 


ds'- 


p(s') 


dxo 


TT^ 


s — ie 2 — t(s', Xg) 


(127) 


The function /(s, A^,m^) is again just the triangle graph of hgure 18 (up to 
conventional constants). 

We are now going to distort the contour C. To do this, we need to 
know about the singularities of as a function of A^. This is 

a rather technical matter, but the most important singularity is easy to 
understand. Looking at hgure 18 along the direction, we can see that 
there is a threshold at = (A + 1)^, which suggests that / has a singularity 
at A^ = (m — 1 )^. 

We can verify the existence of the A^ = (m — 1)^ singularity directly from 
the representation (127). We hrst rewrite the RHS of (127) as 

1 ds' 1 rt+is',m^) 

vr2 X s'- s-ie {[s'- {m-iy][s'- + - t 

(128) 

using (76) - (78). Here t±(s',m^) are the phase space limits in t for a given 
s' and 

t±{s', m^) = (3 + — s')/2 ± 2p(s', nX)q{s'). (129) 

The rightmost integral in (128) is simply 


/ A^ — f_(s',m^)A 
\ A^ — t+(s', m^) / 


(130) 


which has singularities in s' when A^ = t±{s', m^), which is just the boundary 
of the Dalitz plot in the s' — A^ variables. So the singularities in s' are at 
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s'_(A^, rn?) 



.(A^, m^) 


Figure 35: The motion of the singularities sA(A^,m^) for (m + 1)^ < < 

(m — 1)^. 

s'^(A^,m^), which are the intersections of a hxed A^ line with the boundaries 
of the s' — A^ plot. The location of these singularities is what we need to 
understand. 

They can be visualised from figure 14, if we mentally replace s by A^ and 
t by s'. In particular, in the crossed region A^ > (m + 1)^, the intersections 
s'j_(A^, m^) are both negative, and do not interfere with the integration region 
s' > 4 in (128). Thus for large A^ > (m + 1)^ we run into no A^ singularities 
of /. 

Now consider reducing A^. At A^ = (m + 1)^, the points s'j_(A^,m^) 
coincide at the point s' = 1 — m, and then for (m — 1)^ < A^ < (m + 1)^ 
they become complex, one with a positive imaginary part and one with a 
negative imaginary part, as shown in figure 35. The two intersections meet 
again when A^ = (m — 1)^, at the point s' = m + 1, which for m > 3 is 
beyond the start of the s' integration in /. One of s'j_(A^,m^) approaches 
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Figure 36: The original contour C in (124), and the distorted contour. 

s' = m + 1 from just above the contour, the other from just below. This 
means that the contour is “pinched”, which is why there is a singularity of 
/ at = (m — 1)^. 

Note, now, that for m > 3 the singularity at = (m — 1)^ lies to the 
right of the point A^ = 4, so we need to know whether it lies above or below 
the A^ > 4 cut of 4)(A^,m^). The answer is that the physical limit for this 
decay process is taken in the sense of + ie. We therefore position the 
branch point at A^ = (m — 1)^ above the A^ > 4 cut, as shown in figure 36. 

We now distort the contour C so as to wrap around the A^ < (m — 1)^ 
cut as shown in hgure 36. We end up, at this stage, with (126) replaced by 

^(m—1)^ 

2M{s) / dA24)(A^m2)Ai(A^m^s) (131) 

J — CO 

where Ai(A^,m^,s) is l/2i times the discontinuity of across the 

A^ = (m — 1)^ cut. 

The discontinuity of / across the A^ = (m — 1)^ cut can be calculated 
in various ways. Just as we guessed the existence on the A^ = (m — 1)^ 
singularity from inspection of hgure 18, we can guess that the discontinuity 
across the associated cut will be found by cutting the graph as in hgure 37. 
The result will then be proportional to the product of (a) the phase space 
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Figure 37: Graphical representation to get the discontinuity of 
across the = (m — 1)^ cut. 


factor for the intermediate on-shell state of one particle of unit mass and a 
second particle of mass with squared c.m.s. energy and (b) the 

s-wave projection, in the m^-channel c.m.s., of the one particle exchange 
process shown in hgure 38. In this exchange, the momenta are such that 

pI=pI = 1, pI = pI = s, and (pi + (132) 


This is basically correct: from standard techniques of Feynman graph analysis 
[35] [38], the required discontinuity is calculated to be 2iAi(s, where 


Ai(A^, s) = 


In 


'r-Vu^ 


(133) 


7rk{s,m^) \Rj^y/u 

where k{s,m‘^) is as in (77), 

R{s, A^, m^) = —m^ + rri^{s + A^) -|- (A^ — l)(s — 1) (134) 


and 

= k{X^ ,m^)k{s,rn?). (135) 

(A technical detail in parenthesis: the physical region for an external kine¬ 
matic variable like s is in the sense of s -|- ie, but for an internal variable like 
A^ it is in the sense A^ — ie. This implies that the relevant discontinuity in 
A^ will actually be the difference “below the cut - above the cut” [38]. It is 
just this discontinuity that we need in figure 36.) 
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Figure 38: One-particle exchange graph entering into the A^-discontinuity of 
/(s, A2,m2). 

On the other hand, the s-wave projection in the m^-c.m.s. of the exchange 
process of hgure 38 is 


q;(A , m , s) = - dy 


2 7-1 


(136) 


where y is the cosine of the angle between and p^ in the m^-c.m.s. We 
obtain 


a(A^, s) = 


m 


R-VU 

k{s, m^)k{X‘^, w?) """ I i? + ^/u I ' 


In 


(137) 


Further, the two-particle phase space associated with the pi — p 2 state is 
proportional to 

a(A,m) =-5-. (138) 




So we see that, as expected, the function Ai is proportional to the product 
cr(A^,m^)Q;(A^,m^, s). 

Ai will hgure prominently in what follows. It clearly “belongs” in the 
(i.e. 3-body) channel. Up to a kinematic factor, it is the projection of a 
single-particle exchange graph, with the unusual feature that its pole occurs 
in the physical region. For this reason it is often called a “real particle 
exchange (RPE)” process - meaning that the propagator in (136) can vanish 
in the physical region.This singularity shows up in (137), which is singular 
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when = U. This can be written as 4m^r(s, A^, m^) = 0, where T is our 
old friend the Kibble cubic of (80). So Ai has logarithmic singularities on 
the boundary of the s — decay region V. Inside this region, Ai develops 
an imaginary part of ivr. As a result, <h will carry a phase which is additional 
to that of the two-body amplitude M, which is supplied in the isobar model. 
This additional phase is a direct consequence of RPE processes in the three- 
body problem. 


5.6 The single variable representation for $ (2) 

The treatment of the A^ = (m — 1)^ singularity of /(s, A^, vR) was relatively 
simple, but there are other singularities of / at A^ = 0, in this equal mass 
case. These were studied in [35], [38] and [39]. However, in more general 
cases involving non-zero angular momentum states, and particles with spin, 
a simple Feynman graph interpretation is not available. Instead, as men¬ 
tioned earlier, Pasquier and Pasquier [36] showed how the single variable 
representation (SVR) can be derived directly by manipulating the double 
integral in (97). 

In the Pasquier method, one begins by rewriting (97) as 


<h(s, rrX) 


X 


1 1 / 

M(s)+2M(s)- X 

71 s’ -s 

1 

{[s' — {m — l)^][s' — {m + l)2]}i/2 


rt+{s',m?) 
’ t- (s',m^) 


$(t)dt(139) 


where t±(s',m^) are the phase space limits in t for a given s' and m^: 


t±(s', m^) = (3 -I- — s')/2 ± 2p(s', m'^)q(s'). (140) 


The method proceeds by casting the double integral in (139) into the form 
of two contour integrals, and then inverting the order of the integrations via 
a series of contour deformations. The result is that (139) can be transformed 
into the single variable integral equation 


<h(s,m^) = M{s) + 2M{s){ A(A^,m^,s)<h(A^,m^)dA^-|- 


-I- / L(A^, s)<h(A^, m^)dA^}, 


( 141 ) 
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Figure 39: Diagrammatic representation of the integral equation for <F(s, m?) 
omitting the < 0 contributions. 


where 


A(A^, s) 


1 rs+{X'^,m^) 

7r/_(A2,m2) s'-s{[s'-(m-l)2][s'-(m + l)2]}V2 


and 


L(A^, s) = 


1 ds' 


TT 


s{[s' 


[m 


1)2] [s' — (m + l)2]}i/2 

(143) 

Here s±(A^,m2) are the intersections of the phase space boundary curve 
r(s, = 0 with lines of hxed A^. They may be visualized from hgure 

14, redrawn in the variables s, A^ instead of s, t. In the expression for L, 
s_(A^,m^) is the intersection with the boundary of region III. The integral 
in (142) can be evaluated analytically [38] to show that A is in fact the same 
as Ai. The expression for L can also be evaluated [40] in terms of similar 
functions, but we do not give the formulae here. Actually, in sections 6 
and 7 we shall give reasons for omitting the A^ < 0 contribution in (141). 
In this case, the integral equation (141) can be conveniently represented in 
diagrammatic form as in hgure 39. 

The Pasquier inversion was applied to the three-pion system by Pasquier 
and Pasquier [41], and to hnal states of the type KttN (unequal masses, 
non-zero spin) by Brehm [42] and by Aitchison and Brehm [43], [44], [45]. 

All the foregoing can be straightforwardly extended to various more com¬ 
plicated situations. Consider, for example, a model in which we have two 
pairs of hnal state particles interacting so as to form (diherent) isobars, but 
not the third pair. Then we write 


F(s,f, m^) = Mi(s)(/)i(s,m^) M 2 (t) 02 (C 


(144) 
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where the discontinuity across the normal threshold in 0i is 


1 

disc = 2ip(s)-y da;iM2(t)02(t,(145) 

and a similar equation for disc <p 2 - As before, we derive a single variable 
representation for 0i and 02 having the forms 


0i(s,m^) 

— C'i(m^) -|- J 

[ Ki2{t',m‘^,s)(l)2{t',m^)dt' 

(146) 

(t)2{t,m^) 

= (7*2(m^) -|- J 

[ 772i(s',m^,f)0i(s',m^)ds' 

(147) 


where the inhomogeneous terms have now been chosen to reproduce the 
unmodihed isobar terms in the two channels. 

However the equations (146) and (147) are not convenient for practical 
applications, because we would like to be able to calculate the corrections 
from a knowledge of the two-body interactions alone, independent of the 
“htting functions” Ci and C 2 - This is easy to arrange. If we iterate equations 
(146) and (147) we hnd that the equation for F of (144) can be rewritten as 

F(s,f,m^) = [Mi(s)0ii(s,m^)M2(t)02i(t,m^)]Ci(m^) 

[Mi(s)0i2(s,m^) + M2{t)(t)22{t,m^)]C2{rn^) 

= + C2{rri^)'^2{s,t,w?). (148) 

The quantities Tj describe rescattering starting in pair i, and ending in either 
pair i or in pair j. The describe corrections to be applied when an isobar 
is hrst produced in pair j and rescatters dually to pair i. These functions 
satisfy equations of the following form: 


0 ii(s,m^) 


Ki 2 {t', m^, s) 02 i(t', m^)dt' 

(149) 

021 (t,m^) 


{s', t) 0 ii(s', m^)ds' 

(150) 

0 i 2 (s,m^) 

- 1 ^- 

( 0 ,m^,s) 022 (t',m^)d 0 

(151) 

022(0 m^) 


K 2 i{s', t) 0 i 2 (s', m^)ds'. 

(152) 


56 



These coupled integral equations only depend on the two-body amplitudes 
Mj, and can be solved once and for all, leaving the Ciijn?) to be fitted to 
data. 

It is clear that before any of this can be applied to experimental data, we 
must address the complications of isospin and angular momentum. Both of 
these were introduced in a general way into this formalism by Pasquier and 
Pasquier [41]. We now provide a brief introduction to these complications, 
and describe some calculations in physical systems. 

6 Some Practical Examples 

We do not want to get too bogged down in the minutiae of 3-particle helicity 
states - which are contained in the references to be cited. We’ll just give the 
general idea in the case of three spinless particles of unit mass. 

The first step [40] is to generalise the expansion (86) by writing 

= ^(2J + 1){ 

j 

AiU 

A2Z2 

+ Z ®£o(f*3)(24 + (153) 

As^s 

where A* are helicity labels (Cook and Lee [46], Branson et al. [47], Berman 
and Jacob [48]), dij is the angle between the momenta of particles i and j in 
the j — k c.m.s., and hlj specifies the orientation of pi in the 3-body c.m.s. 
Also, li is the pair partial wave and J is the total angnlar momentnm. 

We note that three different complete sets of states are employed in (153), 
so onr basis is overcomplete. However, in each snbenergy channel, only a 
hnite number of partial waves will be retained, and one may regard the low 
partial waves in channels t and u as representing in some average sense the 
omitted high partial waves in channel s. In any case, the reqnirement of 
two-body (subenergy) nnitarity will be imposed, and we shall see to what 
extent three-body nnitarity can be satished also. 
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The second step [40] is to write down the subenergy unitarity relation, 
which takes the form 


+ 2ip{s)M[^*{s)J2i [ C'Ai«iA2Z2*^2'^'''(^,"^^)dcOS012 

A 2 I 2 ^ 

+ similar contribution from $3 (154) 

where is proportional to the Wick [49] recoupling coefficient. If only 

the hrst term on the RHS of (154) were present, we would have the solution 

(s, m^) = (s, (s) (155) 

where G does not have the s > 4 cut, but might have various kinematical 
factors (for example, centrifugal barrier factors as discussed by von Hippel 
and Quigg [50]). This is then the traditional isobar model. The rescattering 
corrections are incorporated by generalising (155) to 

(156) 


The third step is to write a dispersion relation for 0*, and transform it 
into the single variable form. As an example, we write down the equation for 
the = 1“ three-pion (ta) channel (Aitchison and Golding [51] , Aitchison 

|40]): 


(j)^{s,m‘^) = 1 + 2/ dA^Ai^(A^,m^,s)M(A^)0a;(A^,m^) 

J — 00 

+ contributions from A^ < 0 (157) 


where 


Ai^(A^m^s) = 


2k‘^{s, w? 


-(r(s, f, u, m^)Ai(A^, m^, s) — Rk{\^, m^)/2m^) 


(158) 

and T is the Kibble cubic function stu — {w? — 1)^, Ai is as given in (133), 
and R is given in (134). Isospin recoupling is included in (157). 

This equation (including the A^ < 0 pieces) was studied in detail by 
Aitchison and Golding [51]. We parametrised M{s) as 


M(s) 


a + hq^ + cq‘^ + 


2g3 




In 


^/s + y/s - A ^ 

yi- V(s-4)y 


(159) 
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where the log has an imaginary part of —tt for s > 4, and where as usual 
s = 4 + 4g^. Note that the I = 1 partial wave requires the threshold factor 
q^. The parameters a and b are equivalent to the mass and width parameters 
of a B-W amplitude, and were chosen to £t the physical p meson in the first 
instance. We also explored other values of the mass and width parameters, 
as an exercise. The parameter c controls the convergence properties of the 
integral - or in other words the importance of the far left hand region < 0. 
Similar calculations were reported in Pasquier’s thesis [52], but unfortunately 
remain unpublished. 

We found that with nip = (766 — | i 133) MeV it was possible to dynam¬ 
ically generate an u resonance at the physical value, but this required a 
small value of the parameter c, resulting in a strong dependence on A^ < 0 
contributions. We regarded this as unphysical: these left hand contribu¬ 
tions can be thought of as mimicking short-range effects (as opposed to the 
long-range single pion exchange processes associated with the rescatterings), 
which originate in qq dynamics, not 37r dynamics. Nevertheless, it was inter¬ 
esting that the machinery could actually generate a 3-body resonance, and 
the calculated width was satisfactory (but presumably more or less hxed by 
the phase space). 

These considerations led us to try omitting the short range A^ < 0 part, 
so that our equation now reads 

^(m—1)^ 

(j)^{s,m^) = 1 + 2 Ai^(A^,m^,s)M(A^)0^(A^,m^). (160) 

Jo 

This kernel function Ai^ is equal to the appropriate projection of the one-pion 
exchange graph in np —)■ np, up to multiplicative kinematic factors, which 
means that (160) does include all long-range rescatterings. The parameter 
c is now not needed for convergence, and is set to zero, while a and b are 
chosen to fit the p mass and width values. With this truncated equation, we 
expect - and find - much less effect in the mJ channel (no u resonance), but 
pretty much the same result as far as the s-variation is concerned, which is 
dominated by the long-range rescatterings.^ The logarithmic singularity Sb 
mentioned earlier is visible. The deviation in the magnitude of 0^ away from 
unity was generally of the order of 20 - 30 %. A phase of some 20° could be 
generated at s-values in the vicinity of the p resonance. 

^This provides another reason why we may reasonably truncate the A^-integration at 
A^ = 0: the contributions from A^ < 0 are sensibly constant in s, and so may be absorbed 
into the “production vertex” C{mJ). 
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The (oi) wave was also investigated by Pasquier [52] and by 

Parker [53], using the full equations, but the results are only available in 
these authors’ theses. Three channels were included: Tip in both s and d 
waves, and Tre in / = 1 , where e was taken to be a broad low-mass isoscalar 
/ = 0 TT — TT state. Substantial variation was found, as well as signihcant 
rescattering from Tip to Tie. The Svr problem was taken up again by Brehm 
[54] [55], who formulated and solved the integral equations for the = 0“ 
Tip wave, and the l’'' 7 rp and Tie channels. He found substantial w? dependence 
in the case, conhrming the calculations of Pasquier [52] and of Parker [53]. 

The formalism has also been applied to hnal state interactions in tiN —)■ 
TiTiN (Brehm [42], Aitchison and Brehm [43] [44] [45]). The spin of the 
nucleon is a significant technical complication. However, the same steps 
“isobar-type expansion -|- subenergy unitarity -|- analyticity —)■ single variable 
integral equations for the correction functions” can be .followed through. 
The = I"*", \ , I"*", I states were treated. All isobar states likely to be 
important for total energy IT < 1.5 GeV were included, namely Nti isobars 
S'!!, Asi, Pii and P 33 , and tttt isobars in s-wave / = 0, 2, and the p-wave 1 = 1. 
The full integral equations were formulated, but only the first iterations (i.e. 
triangle graph contributions) were calculated, since experience had shown 
that the bulk of the subenergy variation (though not the m? variation) is 
well accounted for by the triangles. 

The main conclusions were as follows. First, none of the corrections 
vary rapidly with the subenergy variable. In addition, the shape of the 
subenergy variation changes very smoothly as W varies. This implies that 
although there may be some observable corrections to the subenergy spectra, 
their presence will not signihcantly distort extracted IT-channel resonance 
behaviour. Thus the non-unitary isobar model was to a large extent justihed, 
at least as the data then stood. That is not to say, however, that with 
vastly more data, and with a focus on interferences on the Dalitz plot, such 
corrections can continue to be neglected. 

Secondly, and in this connection, a number of characteristic subenergy 
variations were found - none very large, to be sure, but possibly signihcant 
nowadays. One such variation exhibited strong curvature at the subenergy 
threshold, the real and imaginary parts crossing over each other. This be¬ 
haviour was found in cases where all orbital angular momenta (both / and 
L) were zero. A simple parametrisation of this pattern is provided by the 
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’scattering length” form 


1 


( 161 ) 


1 — iaqi 

where qi is the magnitude of the pair momentum in their c.m.s. The real and 
imaginary parts of (161) cross at g* = a~^. Typical values of a were in the 
range 0.5 - 1 fm. The zero angular momentum cases were re-examined by 
Brehm[56], who solved coupled integral equations of the type shown in (149) - 
(152) for the relevant amplitudes. The results for the isobar correction factors 
were quite well represented by (complex) scattering length parametrisations. 

A second characteristic variation occurred in which there was curvature 
near the maximum of the kinematically allowed region in the subenergy vari¬ 
able Si- Equivalently, this is the same as peaking in the variable Pi, where 

Pi = {[W - {y/^i + mif][W - {y/s'i-mif]}/2W (162) 

is the magnitude of the momentum of the isobar in pair i, in the 3-body c.m.s., 
and rrii is the mass of the remaining third particle. A simple parametrisation 
of this effect is the form 

-^-r- (163) 

(l + (p,i?)^)2 

We found that R was typically of order 1 fm. Actually, just such a factor 
is frequently introduced into isobar model analyses (along with a threshold 
factor (piR)^), as discussed by von Hippel and Quigg [50]. That such factors 
can have an impact on the subenergy spectrum was noted by Longacre [57]. 

These calculations were done a good many years ago, but were never to my 
knowledge ever combined with a revised isobar-model £t to data. However, 
the equations implementing subenergy unitarity and analyticity are in place, 
and it has recently been stated (Battaglieri et al. [58]) that “with the much 
larger data sets available today, this issue is certainly worth revisiting”. 

At several points in the foregoing the reader will have noticed that, al¬ 
though the initial thrust of the procedure was very much focused on the 
two-body subenergy channels, the end result apparently had relevance to the 
three-body channel as well. The kernel functions Ai are clearly three-body 
in nature, and the rescattering series generated by iterations of the single 
variable integral equations obviously contain three-body intermediate states. 
The question then arises: to what extent do these equations also incorporate 
three-body unitarity? Traditional (Faddeev-type) approaches to three-body 
f.s.i. would of course have a starting point which automatically satisfies 
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three-body unitarity. So it is fair to ask whether the present treatment, 
based on the pair channels rather than the three-body channel, is capturing 
all the relevant physics. In fact, we’ll now see that, rather surprisingly, our 
amplitudes do (or at least can) satisfy an appropriate form of three-body 
unitarity. 


7 Unitarity in the Three-body Channel 

We return to the simple model of (141), retaining only the A = Ai kernel 
(the other terms will not affect the following argument): 

<h(s, m^) = M(s)-I-2M(s) / dA^Ai(A^, s)<h(A^, m^). (164) 

J —oo 

We are now interested in the w? behaviour of F as given by (86). We shall 
take C{w?), the isobar production amplitude, to have no singularities in m?, 
so that whatever three-body structure there is, is in <h(s,m^). 

We hrst verify that <h(s, w?) has a singularity at the three-particle thresh¬ 
old = 9. Note that for > 9 the upper limit of the integration in (164) 
will he to the right of the threshold A^ = 4, where <h(A^,m^) has a branch 
point and associated cut. The physical limit is via the w? F ie prescrip¬ 
tion, and in that case the A^ integration contour will he above the A^ > 4 
cut. If, instead, we give rn? a negative imaginary part, rn? — ie, the A^ 
integration contour will he below the cut, and the result will be different. 
Thus the function must have a branch point at m? = 9, with a discontinu¬ 
ity <h(s+,m^) — <h(s+,m^), which we shall now calculate (we are taking the 
physical limit for s). As usual, this discontinuity will be directly related to 
unitarity in the three-body channel. We shall only aim to give the flavour of 
the analysis: a more careful discussion is contained in Aitchison and Pasquier 
[59], and an even more careful one in Pasquier and Pasquier [36]. 

From (164) it follows that, for tv? > 9, 


<h(s+,m+) — <F(s+,m?_) 


2M{s+){ 


-(m+-l)2 


'r+ 


dA^Ai(A^, s+)<F(A^, m\) 



Al(A^m^s+)$(A^m!)} 


(165) 
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(m+ — 1)' 


A' 



Figure 40: The contours r_ and r+ in (165). 


where F± are the two integration contours lying above and below the > 4 

cut, as shown in hgure 40. We have omitted the + or - labels on the and 

arguments of Ai, since a more careful study shows that the singularities 
of Ai lie on the same (lower) side of both F+ and F_. However, we must use 
<h(A^,m+) on F+, and <F(A?_,m?_) on F_. Thus the RHS of (165) is 


2M(s+) 


“(m—1)^ 


dA^Ai(A^, s+)[<F(A^, m^) — <F(Ai,mi)] 


= 2M(s+){ 




dA^Ai(A^m^s+){[<F(A+,m+) - <F(A+,m?_)] 

+ [4>(A^, m^) - <F(A2_, m2_)]}. (166) 


It follows that 

^(m—1)^ 

disCm2=9<F(s+, m^) = 2M(s+) / dA^Ai(A^, m^, s+)disCm2=9<F(A+, m^) 

J —CO 

+2M(s+) / dA^Ai(A^, m^, s+)discA2=4<F(A^, m?_) (167) 

J —OD 

where discA 2 = 4 <F(A^, m^) means the discontinuity in 4) across the A^ > 4 cut, 
with the prescription — ie. This latter is just the discontinuity given by 
the subenergy unitarity relation, but continued round from m? + ie to m? — ie, 
namely (c.f. (83)) 

discA2=44)(A^, m^) = 2ip(A^)M(A^)F°(A^, m^)6'(A^ —4) (168) 
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where as before = <h(A^, m^)+ <h(t, m^)+ $(«, m^), u = 3 + m^ — 

— t, and is the s-wave projection in the A^-channel c.m.s. of F: 

1 

F°(A^,m^) = - / F{X'^,t,m'^)dxx2. (169) 

2 J-i 


So we arrive at an integral equation for disCm2=9‘h(<s+, m^): 

discm2=9<l>(s+,m^) = 

/ (m— 1 )^ 

dA^Ai(A^, m^, s+) 2 ip(A^)M(A^)F°(A^, m^) 

+2M(s+) / dA^Ai(A^,s+)disCm2=9<h(A^,m^). 

J —oo 


(170) 


The first iteration of this equation is just the inhomogeneous term, which 
we write as 


-(m— 1 )^ 


discm2=9$(s+,m^) = 

dA^2ip(Ai)(T(Ai, mi)F°(A^, m^)\h^^^(Ai, mi, s+) (171) 


where 


>I<0>(A^m^s) = = Z(\\m\s) (172) 

<7{X^^ m^) 


and a was introduced in (138). Now we saw in section 5.4.1 that the quantity 
Ai/a is proportional to the s-wave projection, in the three-body c.m.s., of 
the one-particle exchange process “A^ -|- 1 —)■ 1 -t- s” (see hgure 38), and so 
Z represents just the s-wave projection of hgure 41, which is clearly the hrst 
term in a three particle —)■ three-particle scattering amplitude. Furthermore, 
the factor p(A^)(t(A^, m^) represents the expected phase space factor for the 
effective three-particle phase space (up to conventional constants). We recall 
from (81) that this phase space factor is proportional to 


dXMt 


which, using (76) with s replaced by A^ can be written as 


(173) 


dA= 




2p{X^^m^)q{X^)dx\2 = a{X^^rn?‘)p{}?)dX^ dx\2. 


(174) 
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Figure 41: Single-particle exchange diagram in the 3 —?• 3 amplitude. 

This is just the phase space appearing in (171). The additional 2 in (171) 
arises from the identical channels. 

Carrying out the iteration of (170) to all orders, we hnd the result 


disCm2=9*&(s+>"^^) 



(176) 


where 



d/i^Ai(/i^, s)T(A^ p^). 

(176) 


Equation (175) is a special case of the general discontinuity formula across 
a three-body cut, as given by Hwa [60] and Fleming [61]. These disconti¬ 
nuities were in turn derived from three-body unitarity relations, together 
with subenergy unitarity relations. So it appears that in our approach, a 
combination of two-body unitarity, analyticity, and crossing have generated 
three-body unitarity in the “pair interactions only” approximation, and have 
also generated self-consistently a three particle to three-particle scattering 
amplitude. Naturally our amplitudes contain no three-body forces. Never¬ 
theless, once having made that assumption, it is a viable option to work in 
the two-body channels, rather than in the generally more difficult three-body 
one, without sacrihcing three-body unitarity. 

The amplitude T has a further interesting interpretation. By considering 
the iterative solution of (164), we hnd that <F(s,m^) can be written as 
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where 


s) = Ai(A^,s)+2 / d/i^Ai(/i^,/i^). 

J —oo 

(178) 

Equation (177) shows that the amplitude ip is the resolvent of the integral 
equation for $, where 

ip{X‘^,771^,3) = (t(A^, m^)d'(A^, s)/[M(A^)M(s)]. (179) 

We need to add one important caveat concerning 3-body unitarity in 
this approach. While our amplitudes do satisfy the correct discontinuity 
relations, this only translates into a true 3-body unitarity relation if the 
3 —)■ 3 amplitude is symmetric, as pointed out by Pasquier and Pasquier 
[36] [41] - or at least can be made so by multiplying by a suitable kinematic 
factor. This is the case for the kernel function Ai, but it is not the case 
for the kernels associated with the A^ < 0 part of the integral equations, 
which we have hidden from the reader. One can, indeed, symmetrise all the 
kernels by hand, but then the singularity structure of the amplitude <h will 
be changed. In fact, <h will develop extraneous left-hand cuts in s which are 
unphysical. The most attractive option is to simply truncate the A^ integral 
at A^ = 0, retaining only the Ai kernel - a strategy we already suggested in 
section 6. In this case, of course, one is again introducing new singularities, 
but they lie in the region s > (m-|-1)^, well into the inelastic region [36]. One 
way of seeing that this is likely to be the case is to consider the particular 
rescattering graph of hgure 42. If we want to reconstruct this process exactly 
by dispersing in the variable s, we need to include the discontinuities across 
not only the 2-body cut labelled 1, but also across the 4-body cut labelled 2. 
In the 37r case in the u and Oi channels, satisfactorily symmetric kernels are 
found (Aitchison [40], Pasquier and Pasquier [41]). The same is true for the 
corresponding kernels in the meson -|- meson -|- nucleon case as was shown 
in Aitchison and Brehm [43], where their connection with the elementary 
one-particle exchange processes is also established. 

The upshot of this discussion (see also section 6 and [36]) is that the best 
practical form of the S VR for $ is likely to be the truncated one in which the 
A^ < 0 contributions to the kernel are omitted. We remind the reader that 
in this case the equation has the diagrammatic representation of hgure 39, 
clearly exhibiting the RPE graph in the kernel. 
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Figure 42: Rescattering amplitude. 


Although it is clearly desirable to have amplitudes which do satisfy three- 
body unitarity, in practice this constraint is usually implemented in phe¬ 
nomenological applications as a “quasi two-body” unitarity constraint, ap¬ 
plied to the quasi two-body system consisting of an isobar and a third parti¬ 
cle. Indeed, that was the procedure alluded to in section 3, in final states such 
as 7iK*,pK, etc. In the sorts of systems we are here interested in, namely 
those in which strong isobars are formed, the particle -|- resonance thresh¬ 
old is much more significant physically than the (uncorrelated) three-body 
threshold. But we should understand how this “woolly” quasi two-body uni¬ 
tarity is justified, at least in our model, and how it relates to the three-body 
structure in the <F functions. This will be our last application of the single 
variable representation for $. Such “woolly cut” discontinuity relations were 
discussed by many authors (Nauenberg and Pais [62], Baz [63], Ball, Frazer 
and Nauenberg [64], Zwanziger [65] and Frazer and Hendry [66]). 


8 Particle-Resonance Scattering 

The reduction of the three-body scattering problem to an effective two-body 
one by invoking the isobar model was first introduced by Mandelstam et 
al. [67]. In this approach, we first define what we mean by a “particle -|- 
resonance” amplitude. In the case of <F(s, m^), for example, where <h(s, m?) = 
M{s)(j){s,rri^), we shall take M(s) to have a resonance pole at s = sc on the 
second sheet reached from the upper side of the s > 4 cut, as usual; and we 
also know that 0(s, m^) has the same s > 4 cut but no resonance pole. Let 
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Figure 43: Production amplitude for particle + resonance. 


us denote a value of s on the second sheet of that cut by in Zwanziger’s 
notation. Then M{s^^) may be parametrised near the resonance as 

~ fi'V(sc - (180) 


where sq is complex with a negative imaginary part. We dehne our particle 
resonance amplitude 0(sc, as the residue of the pole in <F(s, m^) at = 
Sc, divided by g‘^\ 

())(sc,m^)= lim ^ ^ <F(s^\m^). (181) 


We represent (f){sc,m?) by hgure 43. In a similar way, from our 3 —)■ 3 
amplitude T(A^,m^,s) (see (176)) we can dehne 


■0(50, Sc) = lim lim 


(sn_sc)(A2n-sc) 


a(A^m^)T(A2",m^s'') 


.11'i 


(182) 

We now turn to the analysis of the singularities of 0(sc,'m^) - in particu¬ 
lar, we want to exhibit the fact that it has a branch point at m? = (y^-|-l)^, 
which is the threshold for the production of the particle -|- resonance state. 
But we must be precise about what sheet this branch point is on. 

We write out our equation for 4)(s,m^) once more, in the form 


<h(s, m^) = M{s) -I- 2M(s) 


y(m-i)2 


dA^Ai(A^,m^,s)M(A^)0(A^,m^) (183) 


where the integration contour in the m? + ie limit lies above the A^ > 4 cut 
of (j) when > 9. Consider now an analytic continuation in starting 
with > 9 and a small positive imaginary part, and then going down into 
the lower half plane, through the real axis. The A^ integration upper limit 
A^ = (m — 1)^ will descend into the second sheet of M(A^) and 0(A^,m^), 
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Figure 44: Calculating the discontinuity across the particle + resonance cut. 

which now become and But M(A^^^) has the pole at 

s = Sc on this sheet, and it follows that there will be a singularity of <F(s, m?) 
in the variable when the end point at A^ = (m — 1)^ hits the pole at 
A^ = Sc- This occurs when w? = + 1)^, the particle + resonance 

threshold (we will see later that it is indeed a square root branch point). We 
can draw an associated cut as in hgure 44. 

The discontinuity across this cut is calculated, as in the case of the 
one across the > 9 cut, by considering the difference between an m?- 
continuation to a point just above the cut (so that the A^ contour T^ leaves 
the pole on its right hand side - see hgure 44), and an m? continuation to a 
point just below the cut (so the contour r_ leaves the pole on its left hand 
side). This difference is 

<F(s+,m^) — <h(s+,m^) = 2M(s+)x 

/ dA^Ai(A^, m^, s+){M(A^^^)0(A^^\ m^) — M(A?_^^)0(A^^^, m^)}(184) 

where now A^^^ (•^-^^) means that the pole in M(A^^^) at A^ = sc is always 
to the right (left) of the oriented contours r+ and r_, and where m\ {m?_) 
means that m? is above (below) the branch cut starting at = (i/sc + 1)^- 
The ± on the A^ argument of 0 is irrelevant since it lacks the pole at A^^^ = sc- 
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The difference in curly brackets is 

The difference M(A^^^) — M[\^_}^) is 27ri times the discontinuity on a full 
circuit around the pole, which is 2nig'^S{X^^^ — sc). Hence (185) becomes 

[<l>(s+,m+) - <l>(s+,m^)] = 2M{s+)27rig‘^(j){sc,m‘i)Ai{sc,m^,s+) 

+2M(s+) / dA^Ai(A^,m^,s+)[<l>(A+“,m+) - $(A+^\m^)]. (186) 

Once again, this is an integral equation for the required discontinuity. 
The first iteration is 


[<l>(s+,m+)-<l>(s+,m?.)](^) = 27ri2M(s+)/0(sc,m^)Ai(sc,m^s+). (187) 

Taking again the residue of both sides at the second sheet pole of M, we 
obtain 


[4>{sc, ml) - 4>{sc, = 27ri 2 g‘^4>{sc, m?_) Ai(sc, sc). (188) 

Referring to (178), this can be written as 

[0(sc,m+)-0(sc,m?.)](^^ = 27ri 2 ^V(sc, ^c) (189) 

where is the hrst iteration of the reduced 3 —)■ 3 amplitude given by 
(179). Iterating the integral equation for our discontinuity then gives 

[0(sc, ml) - (j){sc, m^)] = 27ri 2 /0(sc, m‘^_)^|J{sc, ml, sq). (190) 


Now we know that A^ja is proportional to the s-wave projection of the 
one-particle exchange diagram in the effective 3^3 process “A^ -|-1 —)■ 1 -|- s” 
of hgure 38. So allowing for this kinematical factor, we dehne the reduced 
particle -|- resonance amplitude as 


R(sc,m+,sc) = g^ 


^{sQ,ml,sc) 

a{sQ,ml) 


(191) 


So hnally our discontinuity across the particle -|- resonance branch cut is 


[0(sc, ml) - 0(sc, m‘^_)] = 27ri 2 0(sc, m‘^_)a{sc, m^)i?(sc, sc), (192) 
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Figure 45: The discontinuity across the woolly cut. 


which we can represent as in hgure 45. Note that 

<T(sc,m") = (K - 1)"1K - (193) 

which has the expected branch point at + 1 )^, and is proportional to 

the phase space factor for a two-particle state with masses ^/sc and 1 . 

The extra “2” in (192) arises from the identical t and u channels in this 
toy model. In the case of just a single strong resonance -|- particle state, 
(192) would take exactly the form of (31), with p replaced by a. Thus our 
reduced amplitude 0 (sc,^^) (and of course the similar function <F) satishes 
the expected quasi two-body discontinuity formula, in which the quasi two- 
body —)• quasi two-body scattering amplitude R is self-consistently contained 
in the model, via the repeated rescatterings with one-particle exchange. The 
original isobar production amplitude C{m?), on the other hand, has been 
assumed to have no singularities in m^. 

Similar manipulations show that R{sc, sc) also has the branch point 
at m? = + 1 )^, with discontinuity 

discm 2 =(y^+i )2 R{sc, sc) = 47rii?(sc, sc)o-(sc, m^)R{sc, m^_, sc), 

(194) 

which is the simple (quasi) 2 -body discontinuity relation with p replaced by 
a. We can represent (194) by hgure 46. 

We may regard this aspect of 3-body unitarity as the “long-range” part, 
associated as it is with repeated one-particle exchange between resonances. 
Viewing this dynamically, we might wonder whether such exchanges could 
produce a strong effect in the 3-body channel, perhaps even a resonance. This 
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Figure 46: Discontinuity across the woolly cut of the quasi 2—)-2 particle- 
resonance amplitude R. 


is the same, of course, as asking whether the m^-dependence of the ampli¬ 
tude could ever resemble a resonance. In the case of the truncated 

SVR (i.e. omitting the < 0 part of the integrals), calculations showed 
that this is possible (Aitchison and Golding [51], Pasquier [52], Parker [53]), 
though not for the physically relevant vr — tt amplitudes employed in these Svr 
systems. But in any case, as noted earlier, such physical Svr resonant states 
are generated by qq dynamics, not Stt dynamics. All the same, the extraction 
of the characteristics (e.g. pole positions and residues) of such resonances 
could be substantially affected by m^-dependent rescattering effects. 

It seems natural to place the short-range dynamics leading to resonant 
states in the rn? channel in the production function C{m?). Such states, once 
formed, decay to quasi two-body particle -|- resonance states, and we may 
parametrise such decays by a version of the iP-matrix formalism of section 
3, using the “woolly” phase space function a. We would then end up (in this 
toy model) with an amplitude of the form 

(1 — iKa)~^P [<F(s, m^) + <h(f, m^) + $(«, m^)], (195) 

with a given by (193). Equation (195) represents the formation of an “in¬ 
trinsic” state described by iP-matrix and P-vector poles, decaying as in the 
isobar model but with all pair-wise final state rescatterings included. 
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9 Conclusion 


The main focus of these lectures has been on corrections to the isobar model 
for three-hadron hnal states, so as to incorporate the constraints of two- 
body unitarity and analyticity, primarily, and crossing symmetry somewhat 
indirectly. Initially the formalism appeared to be restricted to imposing these 
constraints on the hnal state two-body channels only. This is in contrast to 
alternative approaches which treat the problem as an essentially three-body 
one, ensuring three-body unitarity as a primary object by the use of model 
Hamiltonians (for example [68]) or via various relativistic scattering theory 
frameworks. 

The “two-body” approach was originally introduced by Khuri and Treiman 
(K-T) [23], in the context of two-particle hnal state interactions in the decay 
K —)■ TTTTTT, with emphasis on the two-body spectra. It is well adapted to 
delivering the functions which correct the isobar model so as to satisfy the 
constraints of two-body unitarity and analyticity, only requiring as input the 
two-body amplitudes (on-shell, but extended as usual in dispersion theory 
beyond the physical region). There is no need for the complicated formalism 
of high spin helds in an ehective Lagrangian, since the usual angular mo¬ 
mentum decomposition of the isobar model is employed. The kinematics is 
relativistic throughout; admittedly, in the case of fermions it causes algebraic 
complications, but these will arise in any approach. It seems fair to say that 
the constraints of two-body unitarity plus analyticity are “minimal”, and 
least model-dependent, in the sense that they should be respected by any 
model. And they turn out to yield a surprising amount of structure. 

Various versions of the K-T approach are available. The one I have 
favoured here is what I call the single variable representation (SVR), where 
the amplitudes obey an integral equation involving only one integration vari¬ 
able. This form may be more convenient for calculations, but its real merit 
is that it is well suited to analysing the model’s properties in the three-body 
channel. Somewhat surprisingly, this simple “two-body unitarity -|- analytic¬ 
ity” model turns out to contain within it a three-body amplitude (consisting 
of pairwise rescatterings), in terms of which three-body unitarity is also sat- 
ished. As originally suggested by Bonnevay [69] [70], this can be traced 
to a form of crossing symmetry, in which the three-body decay amplitude 
“m —)■ 1 -|- 2 -|- 3” is the crossed version (i.e. the analytic continuation in w?) 
of the two-body amplitude “m -|-1 —)■ 2 -|- 3”. The crucial three-body element 
in the SVR is the kernel function Ai (and similar objects in the non-zero 
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spin cases), which is essentially the partial wave projection of the real parti¬ 
cle exchange (RPE) process in the 3 —)■ 3 sector - the quintessential feature of 
three-body dynamics. Arguments were given for retaining only the Ai-type 
kernels in a truncated form of the SVR, represented diagrammatically by 
hgure 39. 

The KT approach is very simple in concept, and delivers isobar model 
correction factors in the form of numerical solutions depending on one suben¬ 
ergy variable (and the three-body mass), given only the two-body amplitudes. 
This technology, developed some 40-50 years ago, can be incorporated into 
fits to three-hadron hnal state physics. At the same time, of course, other 
approaches must continue to be pursued. It will be very interesting to see 
how they compare, in cases where direct comparisons are possible. Perhaps 
we may be able to arrive at a convincingly more accurate version of that 
traditional work-horse, the isobar model. 
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